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Abstract: We compute, in the framework of the fluid/gravity correspondence, the trans- 



> 

On 

fS| port coefficients of a relativistic fluid affected by chiral and gauge-gravitational anomalies, 

including external electromagnetic fields. The computation is performed at first and second 

order in the hydrodynamical expansion. We use a 5-dim holographic model with pure gauge 



^^ and mixed gauge-gravitational Chern-Simons terms in the action. We reproduce at first 

^^ order previous results on the anomaly induced current of a magnetic field and a vortex in 

IL" a relativistic fluid, and compute at second order the anomalous and non anomalous trans- 

port coefficients by using a Weyl covariant formalism. We find a dissipative and anomalous 
correction to the chiral magnetic conductivity due to the time dependence of the magnetic 
^ field. We also find a new contribution from the mixed gauge-gravitational anomaly to the 

shear waves dispersion relation. The role played by the chiral and gravitational anomalies 
in other transport coefficients is discussed. 

Keywords: Gauge-gravity correspondence, Fluid-Gravity correspondence. Anomalies, 
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1. Introduction 

Standard thermodynamics assumes thermodynamical equilibrium, implying that the in- 
tensive parameters (pressure, temperature and chemical potential) are constant along the 
volume of the system. Furthermore it is always possible to find a frame in which the total 
momentum of the system vanishes. In order to study systems in more interesting regimes 
one can allow the thermodynamical parameters to vary in space and time taking the sys- 
tem out of equilibrium. However, we assume local thermodynamical equilibrium which 
means that the variables vary slowly in space and time. This approximation, also called 
hydrodynamical approach, makes sense when the mean free path of the particles is much 
shorter than the characteristic size or length of the system Imfp *C L [1]. 

The modern understanding of hydrodynamics is based on the effective field theory 
formalism. The hydrodynamical systems should obey the (anomalous) conservation laws 
of the spin one currents and the energy-momentum tensor, which are supplemented by 
expressions of the current and the energy-momentum tensor in terms of the quantities in 
the fluid, the so-called constitutive relations. These relations can be written as 

(T^-) = (e+p) U'^U'^ + pg^^'^ + (T^'^)diss & anom , (l-l) 

( J^) = nn'^ + ( J'^)diss & anom • (1-2) 

Here e is the energy density, p the pressure, n the charge density and u^ the local fluid 
velocity. In addition to the equilibrium contributions, there are extra terms in the consti- 
tutive relations which lead to dissipative and anomalous effects. These terms are usually 
computed in the long wavelength approximation, so that they are organized in a derivative 
expansion, also called hydrodynamical expansion. Some examples of dissipative coefficients 
are the shear viscosity r], bulk viscosity C and electric conductivity a (see e.g. [1^3] and 
references therein). 

During the past few years a new set of transport coefficients has been discovered 
as a consequence of chiral anomalies. The axial anomaly of QED is responsible for two 
particularly interesting effects of strong magnetic fields in dense strongly interacting matter. 
At large quark chemical potential fi, chirally restored quark matter gives rise to an axial 
current parallel to the magnetic field [4-6] 

J5 = |^//B, (1.3) 

which may indeed lead to observable effects in strongly magnetized neutron stars and heavy 
ion collisions [7,8]. This phenomena is known as chiral separation effect (CSE). 

In the context of heavy ion collisions it was argued in [9, 10] that the excitation of 
topologically non-trivial gluon field configurations in the early non-equilibrium stages of 
a heavy ion collision might lead to an imbalance in the number of left- and right-handed 
quarks. This situation can be modelled by an axial chemical potential. ^ During the 
collision one expects the generation of magnetic fields that momentarily exceed even those 



^As soon as thermal equilibrium is reached, this imbalance is frozen and it is modelled by a chiral 
chemical potential, at least as long as the electric field is zero. 



found in magnetars. It has been proposed by Kharzeev et al. [9-13] that the analogous 
effect, so-cahed chiral magnetic effect (CME) [14] 

3 = —^^,B, (1.4) 

where J is the electromagnetic current and ^5 the axial chemical potential, could render 
observable event-by-event P and CP violations. Indeed, there is recent experimental evi- 
dence for the CME in the form of charge separation in heavy ion collisions with respect to 
the reaction plane [15,16], and more recently from LHC data [17] (see however [18,19]). 
For lattice studies of this effect, see for example [20,21]. In the context of holography the 
CME was under an intense discussion to confirm its presence at strong coupling [22-26]. 

The fluid/gravity correspondence [27] is a very powerful tool to understand the hy- 
drodynamic regime of quantum field theories with holographic dual. This technique has 
contributed to the understanding of the positivity of the entropy production using tech- 
niques of black hole thermodynamics [28-30]. It is also very useful for the computation 
of transport coefficients. The application of the fluid/gravity correspondence to theories 
including chiral anomalies [31,32] lead to another surprise: it was found that not only a 
magnetic field induces a current but that also a vortex in the fluid leads to an induced 
current, the latter is called chiral vortical effect (CVE). ^ Again it is a consequence of the 
presence of chiral anomalies. It was later realized that the chiral magnetic and vortical 
conductivities are almost completely fixed in the hydrodynamic framework by demanding 
the existence of an entropy current with positive definite divergence [34]. That this cri- 
terion did not fix completely the anomalous transport coefficients was noted in [35], and 
various terms depending on the temperature instead of the chemical potentials were shown 
to be allowed as undetermined integration constants. The contributions from pure gauge 
anomalies is fixed uniquely by this method and provides therefore a non-renormalization 
theorem (see however [36] for a discussion on radiative corrections to the CSE). 

Using a Kubo formula for the chiral vortical conductivity in a system of fermions at 
the weakly coupled regime, a purely temperature dependent contribution was found. This 
contribution was consistent with the integration constants found in [35] and it was shown 
to arise if and only if the system of chiral fermions features a mixed gauge-gravitational 
anomaly [37]. The gravitational anomaly contribution to the chiral vortical effect was also 
established in a strongly coupled AdS/CFT approach and precisely the same result as at 
weak coupling was found [38]. Some evidence of this effect has been found recently also 
from lattice studies [39]. 

Some very recent attempts to establish a non-renormalization theorem for anomalous 
conductivities lead to the fact that the chiral vortical conductivity indeed renormalizes 
due to gluon fiuctuations [40,41]. On the other hand it has been studied in [42,43] the 
ultraviolet cutoff dependence of the anomalous transport coefficients and their holographic 
flow. 



•^A generalization of the model of [31,32] to a Maxwell-Gauss-Bonnet gravity has been done in [33], where 
some corrections of the transport coefficients induced by the Gauss-Bonnet coupling have been computed. 



In [44] the authors claim that the gravitational anomaly produces a Casimir momen- 
tum in the cone formed by the space-time with imaginary time, which breaks the deriva- 
tive counting and it is the responsible that first order transport coefficients being fixed 
by the mixed gravitational anomaly. The gravitational anomaly contribution was con- 
firmed also in a fluid/gravity context [30], in a weakly coupled gas of Weyl fermions in 
arbitrary dimensions [45], and it was found in [46] that the anomalous conductivities can 
be obtained directly from the anomaly polynomial substituting the field strength with the 
chemical potential and the first Pontryagin density by the negative of the temperature 
squared. Recently the anomalous conductivities have also been obtained in effective action 
approaches [47-55] and using group theory techniques [56]. 

Stability and causality issues of the hydrodynamic equations demand the knowledge of 
second order hydrodynamics [57-59]. A classification of the terms contributing to this order 
was presented in [60]. In this work we compute within the fiuid/gravity correspondence 
the transport coefficients at first and second order in the hydrodynamical expansion, using 
an holographic model which includes both gauge and mixed gauge-gravitational anomalies 
and external electromagnetic fields. 

The manuscript is organized as follow. In section 2 we define our holographic model 
and present the renormalized action and the equations of motion. In section 3 we perform 
a formal derivation of the one point functions for a general Lagrangian, either in consistent 
and in covariant form. We review in section 4 the first and second order hydrodynamical 
Weyl covariant formalism, and present our main results. 

In section 5 we explain the method to compute the transport coefficients within the 
fluid/gravity formalism. We present our result of the transport coefficients to first and 
second order in sections 6 and 7 respectively. Finally we conclude with a discussion of our 
results and an outlook towards possible future directions in section 8. The full expressions 
for the sources and transport coefficients at second order are collected in the appendices. 



2. Holographic Model 



The model we will use here was presented in [38]. We will fix first our conventions. We 
choose the five dimensional metric to be of signature (—,+,+,+,+). The epsilon tensor 
has to be distinguished from the epsilon symbol, the latter being defined by e{rtxyz) = +1 
whereas the former is defined by e abode = y/--ge{ABCDE). Five dimensional indices 
are denoted with upper case latin letters. We define an outward pointing normal vector 
UA oc g -^B to the holographic boundary of an asymptotically AdS space with unit norm 



.A 



UATi = 1. The action is given by 



R+ 12 -\fmnF'''' 



^^MNPQR^^^ (^FjvpFqb + XR^ ENPR"" AQRj] + SgH + ScSK , (2.1) 

Sgh = ^ [ d'^xV^K, (2.2) 

ScsK = -^ / d^xV^XnMe'"'''''^''ANKpLDQKJi, (2.3) 

where 5^,^ is the usual Gibbons-Hawking boundary term and Da is the induced covariant 
derivative on the four dimensional cut-off surface. The second boundary term Sqsk was 
motivated in [38] (see also [43]). Notice that the action is diffeomorphism invariant, the 
Chern Simons terms are well formed volume forms and as such they are diffeomorphism 
invariant. They do depend however explicitly on the gauge connection Am- Under gauge 
transformations 5 Am = '^m^ they are therefore invariant only up to a boundary term. This 
model needs a counterterm in order to make the on-shell boundary action well defined '^ 



1 



Srt = — T^^ — ^ I d'^xV—h 



6 + 3P - ( PlfP;: - P2 _ ^F^,F^- ) loge 



where 

6 ' '^ 2 



Bit, - P6i: 



(2.4) 



(2.5) 



Quantities with hat {F, R, . . .) refer to their induced four dimensional objects at the cut-off 
surface, which is located at the radius r ~ 1/e. So taking the limit e — )• 0, one takes the 
surface to the AdS boundary. 

The bulk equations of motion Emn = and M = are 

Gmn + {^F^ - q\ qmn - \fmlFn ^ - 2AeiPQR(Af Vb (F^^i?^ jv) '^'^) = , (2.6) 

y^F^D ^ ^DNPQR ^^F^pFqr + \R^ BNPR'' AQr) = , (2.7) 

and they are gauge and diffeomorphism covariant. 

3. One point functions and Ward identities 

After an ADM decomposition it is possible realize that the action (2.1) is third order in r 
derivatives (see Ref. [38] ) , so in order to get the correct one point functions we have to take 
into account this fact and the assumption that the bulk space is asymptotically anti-de 
Sitter. Asymptotically AdS is enough to get a well defined boundary value problem just 
in terms of the field boundary theory sources. Let us analyze now what this implies for a 
general Lagrange density. 



^It has been proved in [38] that the gravitational Chern-Simons term does not introduce new divergences 
in the system if the space is asymptoticaly AdS. 



3.1 The holographic dictionary with higher derivatives 

Let us assume a general renormalized Lagrangian for an arbitrary set of fields that we will 
call (p. After the four dimensional ADM decomposition one has, 



S 



d^x dr C{(i), (p, Dfj_(f), D^4), i 



where dot indicates derivative with respect to the radial coordinate. A general variation 
of the action leads now to 



SS 



d xdr 



dC,^ dC 



' + 



dC 



AD, 



+ 



dC 



-5{D, 



+ 



dC 



Through a series of partial integrations we can bring this into the following form, 



(3.1) 



5S 



d^xdr^.OM., 



+ 



d^x 



d. 



dC 



D„ 



dC 



d{D, 




The bulk terms are the equations of motion. For a generic boundary, the form of the 
variation shows that Dirichlet boundary conditions can not be imposed. Vanishing of the 
action rather imposes a relation between Scj) and 5(j). 

If we have applications of holography in mind, there is however another way of dealing 
with the boundary term. We suppose now that we are working in an asymptotically anti-de 
Sitter space. The field (j) has therefore a boundary expansion 

(j) = e(^-^>0(°) + subleading , 

here A is the dimension (conformal weight) of the operator that is sourced by (f)^^'. Since 
this is a generic property of holography in asymptotically AdS spaces, we can relate the 
derivative of the variation to the variation itself. 



(A - 4)e 



(A-4)rr,(0) 



subleading . 



Using this and the fact that the one point function of the consistent operator O^ is defined 
as the variation of the on-shell action with respect to the source (f)^^> , we find 



V-/i(o)a 



hm 6(^-4)'- 



— -D (—^^\- — (—\ (A-4)( 



dC 



(3.3) 



Without loss of generality we can evaluate this in Gaussian normal coordinates where the 
metric takes the form ds"^ = dr^ + h^^i^dx'^dx'^, and in the gauge Ar = 0. The gauge 
variation of the action depends only on the intrinsic four dimensional curvature of the 
boundary. From this we can compute the "bare" consistent U(l) current and energy- 
momentum tensor, and the result is 



(c) 



(c) 



-h 



V-/i(o) 






Kt^'^ - Kj'"' + AXe^t^'^^P (^F^pB!;} + Ds{A^R 



Su) 



13 P, 



(3.4) 
. (3.5) 



Now taking the divergence of these expressions and using the equations of motion, we 
get the anomalous charge conservation and the energy-momentum conservation relations 
respectively, 

^M-^f ) = "I^^'^''" (I^M-^PA + Ai?° p,,rP «,a) , (3.6) 

D,Tl^:^ = - J(e)^F^^ + A'^D.J^^^^ . (3.7) 

These are precisely the consistent Ward identities for a theory invariant under diffeomor- 
fisms with a mixed gauge gravitational anomaly. A good general reference for anomalies 
is Bertlmann's book [61] where the consistent form of the anomaly for chiral fermions 
transforming under a U{1)l symmetry group is quoted as 

We use this to fix k and A to the anomaly coefficients for a single chiral fermion transforming 
under a U{1)l symmetry, therefore 

K 1 A 1 , , 

(3.9) 



487rG 967r2 ' levrG 7687r2 ' 

3.2 Covariant form of the current and energy-momentum tensor 

We have computed the currents as the derivative of the field theory quantum action, and 
the anomaly is therefore in the forin of the consistent anomaly. Since we are dealing only 
with a single U{1) symmetry, the (gauge) anomaly is automatically expressed in terms of 
the field strength. However it is always possible to add a Chern-Simons current and to 
redefine the charge current J^ ^ J^ + ce^'^^^AyFpx, and the energy-momentum tensor 
T^^'^ — )• T^^^ + c'e^^^P^Di^ iyAaR^^' px). These redefined quantities can not be expressed as 
the variation of a local functional of the fields with respect to the gauge and metric fields 
respectively. In particular the so-called covariant form of the anomaly difi^ers precisely in 
such a redefinition of the current. ^ 

Adding such a terms to the consistent current and energy-momentum tensor (3.4)- 
(3.5), we can write the covariant expressions for these quantities which are the ones we will 
use to construct the hydrodynamical constitutive relations in the fiuid/gravity approach, 

lQ^GJ^' = f^ F^^L, (3.10) 



MGT^" 



K^" - Kh^"" + 2Ae(^°^^F„^i?^) 



4. Constitutive relations, derivative expansion and Weyl covariance 



(3.11) 



Some notions on conformal/Weyl covariant formalism are needed to construct the consti- 
tutive relations up to second order (for a detailed explanation see [29]). A conformal fiuid 



*Note that the approaches used in [34, 35] and in subsequent works, typicaUy make use of the covariant 
form of the anomaly. 



has to be invariant under the change 



e-'^^^^g,. , 



(4.1) 



where (pix) is an arbitrary function. We will say that a tensor is Weyl convariant with 
weight w if it transforms as 

Qf.::. ^ e-<^(-')Q^^;;; . (4.2) 

The consequences of conformal symmetry on hydrodynamics is that the energy momen- 
tum tensor and (non)-conserved currents have to be covariant under Weyl transformations 
and the energy momentum has to be traceless modulo contributions from Weyl anomaly. 
To construct Weyl covariant quantities it is necessary to introduce the Weyl connection 

1 



wAi, 



U^D^Uy 



and the Weyl covariant derivative 



'DxQi:. 



DxQ'i::, - wAxQii: 

+ [gXaA^" - 5lAa - 

- [gxuA'^ - 5lAu - 



-D^v!" 



+ 



(4.3) 



(4.4) 



We show in table 1 the Weyl weights of some of the hydrodynamical variables. It is 
possible to reduce in a systematic way the number of independent sources contributing to 
the constitutive relations by imposing Weyl covariance and the hydrodynamical equations 
of motion (Ward idetities). A classification in the so called Landau frame of all the possible 
terms that can appear in the energy- momentum tensor and U(l) current has been done 
up to second order in [31,32,60]. The Ward identities in four dimensions in presence of 
quantum anomalies are shown in (3.6) and (3.7). The curvature part has been usually 
neglected in the literature as it is fourth order in derivatives and the expansion is usually 
done up to second order. But it was shown in [37,38] that the gravitational anomaly indeed 
fixes part of the transport coefficients at first order. Actually in [44] it was understood 
why the derivative expansion breaks down in presence of the gravitational anomaly. 



Field 


weight 


^i, T, u^' 

g^iu 

P 
n, E^', B^ 


1 

-2 
4 
3 



Table 1: Weyl weights for the chemical potential, temperature, fluid velocity, metric, pressure, 
charge density, electric field and magnetic field. 



With these ingredients we can write down the constitutive relations in the Landau 
frame 






piAu^u-^ + vn + rf,'; + r^';;,, + r^:. + r,^^ 



(2) 



(2)ano 



nu^ + u'r^. + V 



+ <o^ + z^: 



(l)ano ^ ""(2) ^ '^(2)ano ' 



(4.5) 
(4.6) 



where we have spht the expressions in the equihbrium, first order and second order (anoma- 
lous + non anomalous) parts. Weyl invariance implies the equation of state e = 3p and the 
vanishing of the bulk viscosity C = 0. ^ The subindex in parenthesis indicates the order in 
the derivative expansion. The ambiguity in the definition of temperature, chemical poten- 
tial and fluid velocity which appears when the system is slightly out of equilibrium is fixed 
by using the Landau frame, in which it is demanded that u^^T(^n)fiu = ^^^ ''^^^{n)fi = 0- Up 
to first order, the most general contributions to the conformal energy momentum tensor 
and U{1) current are 

^S = -2W^^ -ri)ano = 0' (4-7) 

i^^,, = -a {TP^^-V^f, - En , <na„o = ^VOJ'' + ^bB^ > (4-8) 



where we have defined the Weyl invariant quantity fl = /u/T. In these expressions ij, a, ^y 
and S,B are the shear viscosity, electrical conductivity, chiral vortical and chiral magnetic 
conductivities respectively, while a^'^, lo^, E^ and B^ are the shear tensor, vorticity, electric 
and magnetic field respectively, defined as 

^/i!^ = 2 (^M^*^ + ^i^'^fj.) ) (4-9) 

(^fiu = ^^^lUu - VyU^ , (4.10) 

u^ = h^'^^^u^u^p, (4.11) 

E^' = F'"'uu, (4.12) 

B^ = ^e^'^f'n^F^i,. (4.13) 

Finally P^'^ = h'^'^ + u'^u'^ is the projector in the space orthogonal to the velocity field. 
From previous definitions one can easily prove that the strength tensor decomposes in the 
following way: 

F^fS = UaEp - UfsEa - e^ppf.uPB^' . (4.14) 

The second order contributions to the constitutive relations are 

a=l a=l 

a=10 a=5 

a=l a=l 



■'As it has been discussed in [60] the anomalous terms in the constitutive relations are those whose 
transport coefBcients have {C,P) = (±1,-1), and they correspond to the ones containing odd powers in 
the anomaly coefficients k and A. See this reference for a systematic classification of anomalous and non 
anomalous terms. 



with the second order tensors defined as 

(4.17) 



f(i)^lu ^ p(/^^i'> ^ f{5)^lu = 5(Mp^)/2 , 7-(6)m;^ = ^(^S^) , (4.18) 



where C'^"'''^ is the conformal Weyl curvature tensor, H'^'^ ap = \ [PaP^ + P^PJ^ - ^P^^Pafi 
is a transverse traceless projector and we use the notation X^f^'^' = U^^'^ ^^X"^. The second 
order vectors are 



' i^a J 



j(io)m = e^'^'^^UyV^BB , 






(4.19) 



(4.20) 



This classification of independent possible terms has been previously presented in [60] . We 
consider in this work a flat background metric, so we are neglecting those terms proportional 
to the conformal Weyl tensor. 

Here we write for completeness the first order transport coefficients for the anomalous 
holographic plasma computed in the literature. 



r 



3 



■' " IGttG ' '^ " WG^ ' ^^-^^^ 

^^ ^ V3q{m + 3rX)K ^ ^^rqT^X ^^ ^ _^^ ^ 27r (2g^ - rj) T^X 
SvrGmr^ Gm ' ATiGm Gmr\ 

where ttt., q and r+ are the mass, charge and radius of the outer horizon of the black 
hole, which we will define in more detail in section 5. We write also here our main results 
corresponding to the second order transport coefficients that are completely new. We 



10 



express some of them as an expansion in the parameter n assuming that /i <C 1, 

As = -^(^-^(141og2-l),^ + 0(,^)), 
A9 = ;l-(y + ^log2 + ^(7-81og2)/22 + 0(/2^; 



(4.23) 
(4.24) 
(4.25) 



^12 - o 






256 

+ ^^kA(5 - 12 log 2) + \^{29 + 60 log 2) + Oijj?^ 

Idtt 



(4.26) 



A5 
Ae 

A7 



A4 



1 



2r? 2ttpG 
Arjfl 



-Ai 



K^ + A(7rr)' 



2^2 v'^(l-21og2) + 2A(l + 21og2) + C'(/i2)) , 
3 ^ -Klog2 + A(5 + 21og2) + 0(/22)) , 



Srjfl 



32r/ 



r 



2 ^' 



As — ^ 



+ 



64(^G)2 
riA 



87r2G2 



5 
18"^' 



(4.27) 
(4.28) 
(4.29) 
(4.30) 
(4.31) 



log2 + -^/22(5 + 61og2) + -^/i4(199 - 7261og2) + 0(/z6' 



1 
2167r4' 



1 np.^ k'^jP log 2 Jl'^\og2 0/ 1456/^2 16/^2 log 2 

?6 = - -. „ ^ — o^ H 5-TT — TT^ ^ + .„ ,^ + A 



+k\ 



levrG 1927r3G tt^G UttG AStt^G \1^t:^G tt^G 

4 130//2 8 log 2 20/2^ log 2^ 



^7 
6 



(4.32) 

(4.33) 

(4.34) 

7i / 7 2, , , 4KA(14-271og2) 4A2(31 - 601og2)\ ^^, o, 
Vt:: ( TTTTT - 2^2(2 - 3 log 2) + ^ ^^ ^ ^-^ + 0{p?) , 



Gtt 9G^3 + Qt^ ^Qt^3 



+ 0(/i'), 



' ^' + ^°^V:^Al-81og2))+0(;2^), 



levrG V 2 127r2 

JL(J_ 
7r3GV192 

1 / 11 



TT^Gr+ \ 192 



+ 2k^ log 2 - 2kA(1 + 2 log 2) + 39A^ fl + ©(/i^) , 



1/^1, //^ 2/i2 2, ^ , ol__o^ , 2/u' 



^A' (9w' - ,1' fll^ + 16 log 2") ) U 0(;i*) , 



TT^ 



15 



(4.35) 
(4.36) 

(4.37) 

(4.38) 
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fj- 



1447r4G 



3K(247r2 log 2 + //^(SS - 68 log 2)) - 8A(187r2(l + log 2) + /i^(5 - 51 log 2)) + 0(/i^) , 

(4.39) 



6 

Ie 



1 ^j21os2+55!P^^i^V4AA+2te2 f^^P + ^'-^^i 



8G7r2 



vr^ 



67r2 



0(/i^ 



^5 



a 



8/2 

7rr+ 



(k log 2 - 2A(1 + 2 log 2)) + 0(/2^) . 



As we shall see most of the coefficients computed receive A— corrections. Indeed not only 
the anomalous transport coefficients are sensitive to the presence of the anomaly, but also 
the non-anomalous ones get corrected as well. The latter had been computed in the past 
without including the mixed gravitational anomaly, so these corrections were neglected. 

5. Fluid/Gravity Computation 

The system of bulk equations of motion (2.6) and (2.7) admits an AdS Reissner-Nordstrom 
black-brane solution of the form 



ds2 = _^2^(^)rfi2 ^ dr 



r'ifir) 



+ dx^dx^ , 



A = <t>{r)dt, 



(5.1) 
(5.2) 



with /(r) = 1 — m/r'^ + q^ /r^ and (l){r) = — vSg/r^ . The real and positive zeros of /(r) 
are 



-+ = ?fl + AA + ^/^' 




-1 + 



\ 



1 + a/1 + 



2 r,2 



(5.3) 



(5.4) 



Stt- 



/i- 



where r+ is the outer horizon and r_ the inner one. The mass of the black hole can be 
written in terms of hydrodynamical variables as 

. 3 / 



m= ""—^ I 1 + ^/1 + ^/22 



-1 + 3a/ 1 + 7^^/22 



24 \ ' ' V ' ' 37r2^ / \ ' ' "V ' ' 37r' 
The boosted version of this blackhole in Eddington-Finkelstein coordinates looks like 
ds^ = -r'^j{r)u^Uydx^dx'' + r'^P^ydx^dx'' - 2uy,dx^dr , 



(5.5) 



A 



{r)Ufj_dx^^ , 



(5.6) 
(5.7) 



with the normalization condition u^«'^ = —1. (5.6) and (5.7) is a solution of the equations 
of motion as long as m, q and u^ are independent of the space-time coordinates x^. The 



(4.40) 

(4.41) 
(4.42) 
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fluid/gravity approach tells us that we have to promote all the parameters to slow varying 
functions of the space time coordinates, and include corrections to the metric in order to 
make it a solution of the equations of motion again. 

5.1 Weyl covariant ansatz 

In order to follow the fluid/gravity techniques [27,31,32,62] we will use a Weyl invariant 
formalism [63] in which the identification of the transport coefficients is direct. We start 
with the ansatz 



ds^ = -2Wx{p)u^dx'' {dr'^ + rAudx") + 



W2{p)ri^y + W^{p)u^Uy + 2 



Wi,{p) 



P^uu 



+ 



W,^,{p) 



dx^'dx'^ 



A 



,W 



+ au{p)Pl^ + r+c{p)uA dxf" 



(5.8) 
(5.9) 



where now r^ is an unknown function of the space-time coordinates, r^{x^), and coincides 
with the radius of the outer horizon of the black hole (5.3) only when the x'^ dependence 
is gone, a^ = a^'(a;'') is a boundary background gauge field satisfying a}t (xq) = 0. 
Notice that r^^i/ is the Minkowski metric, so we will look for metric solutions with flat 
boundary. The r-coordinate has Weyl weight +1, and in consequence r+ has the same 
property. By construction the W functions are Weyl invariant, so that they will depend 
on r only in a Weyl invariant way, i.e, W{r) = W{p) with p = r/r^. W^^^i^r) obeys 
the traceless and transversality conditions W^ {r) = 0, u^W5^u{r) = 0. All these scalars, 
vectors and tensors will be understood in term of a derivative expansion in the transverse 
coordinates, i.e. F{p) = F^^^p)+eF^^\p)+e'^F^^\p)+0{e^) for a generic function F, with 
e a parameter counting the number of boundary space-time derivatives. ^ This solution 
leads to the current and energy momentum tensor after using the AdS/CFT dictionary 
(see (3.10) and (3.11)) 



Jn 



^ lim 
SvrG €^0 



ric(2» 



u^ + rlaf'^ + J^ 



ct 



T, 



1 



jiV 



167rG e^o 



lini [rX{yV2 + W^f^^'\^u^u, + r?^,; 



^^rlwt:} ^^T\wff^P?..u 



5fiu 



'4(T 



(At "i^) + '^^iv 



(5.11) 



where F^^'*"' denotes the coefficient of the term (p~^ — e)" in an expansion around the 
regularized boundary, and e determines the position of the cut-off surface p = 1/e. The 
counterterms in the current J^ and energy momentum tensor T^^ are needed to make the 
expressions finite, and they follow from the counterterm of the action (2.4). They write 



jct 


= 2^08 


e 


{V^E'' - 2uj^B'') 


u,+4'^-l4'^ 


rpCt 


= loge 


" - \{BpBP + EpE^)P^, - \{E^E^ + 




-{jti^'^u^ + j^'K) 





jP + J, 



(10) 



(5.12) 



i?«i?")n^7x, + r(^) + r(3°) 



(5.13) 



^This expansion of F{p) is basically a Taylor expansion around the point Iq . 
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We are considering a flat background metric, and so the divergences appear only 
through terms involving electromagnetic fields, in addition to the cosmological constant 
contribution which was already taken into account in (5.11). 

The functions at zeroth order in the derivative expansion correspond to the boosted 
charged blackhole, i.e. ^ 

c(^p) = -^, (5.14) 

wi'\p) = l = W^'\p), (5.15) 

Wi'\p) = l-f{p), (5.16) 

wi'Jip) = = wt,Up), (5.17) 

aP{p) = 0. (5.18) 

Then the charge current and energy momentum tensor at this order read 



4°^ = £S-. ' T^'J = T^ (4^M^^ + V,.) ■ (5.19) 



From this we obtain the equilibrium pressure and charge density p = jjSj and n = g^. 
For computational reasons it is convenient to define a Weyl invariant charge Q = q/r"^ 
and mass M = m/r'^ = 1 + Q^. In terms of these redefined parameters, the black hole 
temperature and chemical potential read 

r=^(2-Q2)^ ^^^^_^g^ (520) 

ZTT 

We also define the inner horizon in the p-coordinate, p2 = r^/r^ . 

5.2 Einstein-Maxwell equations of motion and Ward identities 

Inserting the ansatz (5.8)-(5.9) into the Einstein-Maxwell system of equations we find a set 
of (2 X 1 + 2 X 3 + 5) differential equations and (2x1 + 3) constraints relating the allowed 
Q{x^), r_|_(x^), u'^{x^) and al, {x^) [32]. ^ We need to solve the equations of motion around 
a certain point Xq that we choose to be Xg = 0. At such point we sit in a frame in which 
u^'{<^) = (1,0,0,0) and aif^(O) = 0. 

The scalar sector is obtained from the rr^ rv and vv components of the Einstein 
equations, and the r and v components of the Maxwell equations. One finds two constraints 

Et^ + r^f{r)Ei^^ = ^ {V^T^ = F,, J")^"-^) , (5.21) 

Mi") + r'^f{r)M';P^ =0 ^ (V^J^' = ciF A F)^"-^) , (5.22) 

which, as indicated, correspond to the energy-momentum and current non-conservation 
relations at order n — 1. ^ The combinations E^-r = 0, Erv + r"^ f{r)Err = and M^. = 



^Following the notation in [31], barred superscripts {n) should not be confused with superscripts (n), 
where the latter refers to the order in the hydrodynamical expansion. 

*1, 3 and 5 denote the 5*0(3) scalars, vectors and tensors in which the fields are decomposed. 

^Notice that there are no curvature terms in (5.22) and (5.21) because we are working with a flat 
boundary. 
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leads respectively to the set of differential equations 

3 



?>dM^\p) - \p-^dp (/dpWtHp)) = §(")(/>) , (5.23) 
dp (/Wf'') + ^p'wf^ - ^Qapc(") + (1 - V)apW^i") - Vv^i"^ = k(")(p) , (5.24) 
dp UdpC^'^A - 2^QdpW^^\p) + ?>^QdpWt^ = C(")(/9) . (5.25) 



At this stage there is still some gauge freedom in the metric. There are three (metric) 
scalar fields (VFi , W2 and Wz) but the Einstein's equations give two differential equations, 
(5.23) and (5.24). We choose the gauge W2{p) = 1 in which the system partially decouples 
and can be solved as 

1 

Wi''''\p) = -\ 1^ dxS(")(x) , (5.26) 

3 J p 

wi^'^'Xp) = ^ - 1 /"' dx f]K(")(x) - ?,x^w[^^'\x) + ^5,c("'^)(rE)) , (5.27) 
P P J p \ v3 / 

c(-.^)(p) = CO ^^ ~^y^ - p dxx-'^ f^ dy (c(")(y) + ^S^"^(^)) • (5.28) 

These solutions have been constructed by requiring Dirichlet boundary conditions at the 
cut-off surface and demanding regularity at the interior of the bulk. The remaining in- 
tegration constants Cq and cq are associated to the freedom of choosing a frame in the 
hydrodynamic set up. 

In a similar way, the vector sector is constructed with the components of the equations 
of motion Eri, Em and Mj. They lead to a constraint equation, 

^i?+^'/(O^i?=0 => {VpTt = F,^rt-'\ (5.29) 

implying the energy conservation equation, and the two dynamical equations 

dp {p'dpW^^ + 2 Vsgaf ) (p)) = jj") (p) , (5.30) 

dp (p'f{p)dpa^\p) + 2V3QdpW^\p)) = A^\p) , (5.31) 



corresponding to Eri = and Adi = respectively. The general solution of this system in 
the Landau frame has been found in [31] with neither background fields at the boundary 
nor gravitational anomaly. It is straightforward to generalize the solution to the case in 
which electromagnetic sources are included. In this case new divergences arise that need 
to be regulated with the cut-off 1/e, and then to be substracted with the corresponding 
counterterms (5.12)-(5.13). Having found in [31] the general solution of the system, the 
contribution to the current coming from the vector sector writes (cf. (5.10)) 

4^'^) = I (-K H/e) - p dxKix)] - ^C(^) - ^Dt\ (5.32) 
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where the integration constants Cu and D^ are determined by fixing the Landau frame 
(^lo-*^ = 0), and demanding regularity at the outer horizon respectively. These constants 
write 

m=0 

D(;^ = -V3Q /' dx M^ - M /' dx ^ + Q2 /' dx ^ . (5.34) 

Ji x^ Ji x'^ Ji x^ 

Finally the tensor equations are the combination Eij — g^ytr {Eki) = 0, which leads 
to the dynamical equation 

dp {p'f{p)dpW^]{p)) = pS;V) • (5.35) 

The solution of this equation that satisfies the Dirichlet boundary and regularity conditions 
writes 

W^f.u'ip) = -J dx ^^^^^^ . (5.36) 

After doing an asymptotic expansion of this solution around the regularized boundary 
surface, one can extract the tensor contribution to the energy momentum tensor (cf. (5.11)), 

,„,(4,.) ^ (-i)-a-p^.(i/6) p.p.. ._„. 

Note that the form of the homogeneous part in the dynamical equations in the scalar, 
vector and tensor sectors is the same at any order in the derivative expansion. Each order n 
is then characterized by the specific form of the sources. In the next two sections we will 
compute the sources, and integrate them according to the formulae presented above to get 
the transport coefficients at first and second order. 

6. First Order Transport CoefRcients 

The technology presented in Sec. 5 can be used to construct the solutions of the system 
at any order in a derivative expansion. As it has been already explained, the solution at 
zeroth order trivially leads to the charged blackhole with constant parameters (5.6)-(5.7). 
In this section we will solve the system up to first order. The transport coefficients at 
this order have been obtained previously in the literature using different methods in field 
theory and holography. In particular, they have been computed within the fluid/gravity 
approach, but not including external electric fields in this formalism, see eg. [30-32,34]. 

6.1 Scalar sector 

In the scalar sector, the first order sources look like 

sW(p) = ]kW(p) = cW(p) = 0. (6.1) 



16 



This 


very 


simple 


situation leads to the solution 










w['- 


'\P) - 


= 0, 












c^'- 


''\p) -- 


^.o(^- 


- e'p') 
P' ' 










wt 


'\P) - 


p'^ 


2Qco (1 
V3 


- e'p') 
P' 



(6.2) 
(6.3) 

(6.4) 

The integration constants cq and Cq can be fixed to zero because they just redefine the 
charge and mass of the black hole respectively. 

6.2 Vector and tensor sector 

The first order sources are given by 

,« . -.^ (f - m) ^ - ^«A (1»^ - ^) .„ , ,0.5, 

(1) _ VSttT n (^^ 9Q2 \ E^ IQ^kQB^ 48kQ2 



r_|_ p3 j.^ ^5 



48A (15Q4 - 16MQV^ + 4MV^) 



pii 



(6.6) 



P« = -6r+pV^, , (6.7) 



where V^ is the Weyl covariant derivative and V^Q = r- ^^ — —VaP- Using equations 

(5.10), (5.11), (5.32), (5.33) and (5.34) it is straightforward to find the first order transport 
coefficients shown at the end of section 4. We write the result again for completeness. 



r 



3 






V3q{m + 3rl) k VSngT^X 3q^K 2^ {2q^ - r%) T^ \ 

^"^ ~ SGmvrr^ + Gm ' ^^ " 4Gm7r + Gmrl '^ > 

Chiral magnetic ^b and vortical ^y conductivities have been computed at first order in 
holography within the Kubo formulae formalism in [38,64-67], including chiral and gauge 
gravitational anomalies. Here we reproduce the same result within the fluid/gravity ap- 
proach. ^° 

Note that to compute the first order transport coefficients one needs only the terms 
ajj, '"" and VF'^'*^) in the near boundary expansion. However, in order to go to the next order 
in the derivative expansion, we need to know the exact solutions, which can be written in 
terms of the sources as 

W^^^ip) = Fi[p\P;V,Q{x) + F2[p\u^{x) + F-M 

W^^^,{p)=F^[p]r+a^,{x), 

a«(p) = Fq[p\P;V,Q{x) + Fj[p\u^{x) + F^ 





(6.10) 




(6.11) 


"^^^^KF.lpf^^^^ 


(6.12) 



r-^ Va 
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The gauge gravitational anomaly contribution to the chiral vortical conductivity was also computed 
recently within an holographic setup in [30,68]. 
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We show in Appendix A the expressions for the F^s functions. F^ writes 



^5[P] 



21og[l + p] _ {1 + P2 + pI) log [p - P2] 2 (1 + pI) log [p + p2] 



-1 + M (i + p2)(l + 2/5i) 



+ 2 + bpl + 2p\ +2 + 5p2 + 2p4 



ArcCot 



-2 - 2pl + 4p4 



P 



vTTpf 



(6.13) 



7. Second Order Transport Coefficients 

The second order coefficients are much more computationally demanding than the first 
order ones. The parameter c*-^'*^^ in (5.10) can always be chosen to be zero, as it just 
redefines the charge and mass of the black hole. On the other hand, because we are 
working in the Landau frame, there is no contribution coming from the scalar sector to the 
energy-momentum tensor and {W2 + VFs)^^'^^ is set to zero. We have checked that this is 
in fact what happens by using the sources for the scalar sector. So, we will focus in this 
section on the vector and tensor contributions. 



7.1 Vector sector 

The second order sources in the vector sector are shown in Appendix B.l. Again using 
these expressions and Eqs. (5.32), (5.33), (5.34) and (5.10), we can extract the second 
order transport coefficients. We show first the new non anomalous coefficients 



?5 = ^5,0 (P2) , 
?6 = '^6,0 (P2) + 

?7 = ^7,0(^2) , 



3 (3 + M2) Q^k 
47rGAf3 



+ i^Kq,k\{P2) + >?i&,Xi{P2) 



+ l^^i?,,K^{P2) + kXC8,kx{P2) + A^C8,A2(/'2) 



(7.1) 

(7.2) 
(7.3) 

(7.4) 



(9 + 12M + 7M^)ttQT^ 
128V3GM^{l + M)rl 

Q (88 + 480Q2M + 169g6) 1 , „ , _ ^ . ^, ^ 9 , w ^^ 

^ ^ 2 + _ {k%^AP2){P2) + kXC9,.x{P2){P2) + X%MP2){P2)) 



512^7rGM4r+ 



r+ 



(4 + 7(5^) 
60 = Q^^Qj^^ + ^^Cioa2(P2) + /«ACio,ka(P2) + A2Cio,a2(P2) 



(7.5) 
(7.6) 



These coefficients had not been computed previously in the literature. The rest of the non 
anomalous coefficients were obtained in the past without the gravitational anomaly. In 
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this work we have found the A— corrected results, which write 



6 



8GM3(M + l)r2 



Q^ + 



1 + 2pl 



log 



1 



Pi 



(7.7) 



(3 + M)(M(3 + M)-6)r2 



GM3(M + l)r2 



6 

^4 



128GM3(M + l)r+ 
3\/3Q3r+ 

„2 

+ r+KA^4,KA(p2) + r+A^^4^A2(/02) • 



U-kGNP ' 
3\/3QV+k2 



(7.8) 
(7.9) 

(7.10) 



27rGM2 
In the anomalous sector, the new coefficients (not computed previously) write 



6 



6 



3\/3Q3 (6 + M) k \/3vrgr2K log 



+ 



2+p| 



16^GM2 (1 + 2/92)2 ■ 2GMr2_ (1 + 2p2) 



,2\3 



+ A|2,a(/02) 



SrlM^G 



Q'^K + 



iTT^T^X 



ii = <4_k(p2) + Af4,A(p2) , 
Is = <5,k(/02) + X^5,X{P2) , 

while the already known coefficient with the new A contribution writes 

3Q2r+K 



6 



+ Ar+^i^A(P2)- 



(7.11) 

(7.12) 

(7.13) 
(7.14) 



(7.15) 



47rGM2 

The ?i,(o,K2 ,jA,a2)(/02) and ^i^[K,x){P2) functions are defined in Appendix C. These coefficients 
enter in the constitutive relation for the current through (4.6) and (4.16). 

7.2 Tensor sector 

The second order sources in the tensor sector are shown in Appendix B.2, and again we 
can extract the transport coefficients at this order after pluging these expressions into 
Eqs. (5.11) and (5.37). Due to the length of the expressions, some of them will be shown 
exactly and the rest are expressed in terms of some functions Aj (^2 ,jA,a2)(P2) and Aj (^ a)(P2) 
which are presented in the Appendix C. Again we split our results in those non anomalous 
coefficients which are new. 



A7 



\/3(-3 + 5M)Qr+ 



647rGM2 
As = r+A8,o(p2) , 

Ag = A9,o(p2) , 



Alo 
An 

A12 
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+ '^^Aio,«;2(/32) + kAAio,ka(P2) + X'^Ainx^{p2) 



96ttG 
sV^Qv+kX 

V3Qr+ 3\/3Q3r+K2 



(7.16) 
(7.17) 
(7.18) 

(7.19) 
(7.20) 



levrG 



+ 



ttGM 



+ KAr+Ai2,KA(/32) + AV+Ai2,a2 {P2) , (7.21) 
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and the rest of the non anomalous ones 

ri /„ M . [3- V4M-3 



"''TEkii'^VSTTt""^ 



3 + V4M-3 



(7.22) 



„2 



'^^ = Sfe • <^-^=« 



„2 



M 



^^ = 8G^V^(rT^^°' 



2 + pI 



pI 



+ 192Q^.A-^^^(^^-^^-^^^, (7.24) 



QM 3QViK^ 18Q^(5 + Q^(9Q^-16))ri.A 

"^^ = "16^ + ^^GiW- + S^^GM + ^ ^+A4,,.(/,2) , 

(7.25) 

A. = ^ -^^^ , (7.26) 

16V3GM2(M + l)r+ ^ ^ 

Ae = riA6,o(/02) • (7.27) 

For the anomalous coefficients we get the new ones 

As = Kr+A5,«(p2) + Ar+A5,A(P2) , (7.29) 

Ae = kA6,k(p2) + AA6,a(p2) , (7.30) 

A,^-^, (7.31, 

3g2 (g2 _ ^\ 

As = "" ZgM^ + Ar+As,,(p2) , (7.32) 

and the rest of the anomalous ones 

~ V^rlQ^K ^ V3Qr+(3r+ + (Q2_4)7rr) ^ 

^' = - At^GM + ^^GM ^ ' (^-^^^ 

A2 = Ar^A2,A(p2) , (7.34) 

97^2 \ 

These coefficients enter in the constitutive relation for the energy-momentum tensor through 
(4.5) and (4.15). 

The transport coefficients Ai, . . . , Ae, Ai, . . . , A3 and ^1 ... ^4, .^1 have been computed 
in the past in [31, 32] without gravitational anomaly. It is interesting to remark that 
Ai, A2, A5, Ag, A7, Ag, Ag, ^i,,^3,,^5 and ^7 do not receive A— corrections, actually these co- 
efficients do not depend on k either. It is also remarkable that A2 and A3 in the presence 
of gravitational anomaly are not vanishing. The rest of the transport coefficients we have 
computed are new. 

7.3 Discussion of second order results 

It would be interesting to compare our results with the predictions done in [60]. Basically 
the authors tried to fix the anomalous second order transport coefficients using a generalized 
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version of the method developed by Son & Surowka [34]. The only issue is that they didn't 
consider the mixed gauge-gravitational anomaly and neglected all the integration constants 
as the previous authors. Nowadays we know that at least at first order these integration 
constants might be related to the anomalous parameter A. The authors presented a set of 
algebraic and differential constraints. The algebraic ones are 

Ai = - i^v - TDb) , (7.36) 

n 

h = - i^B - -^fi) , (7.37) 

n 

Is = —{^v- TDb) , (7.38) 

n 

Is = , (7.39) 

where Db = j^A*^ is the coefficient multiplying the magnetic field in the entropy current 
computed in [34, 35, 60] with only pure gauge anomaly, and R = 2^ is the anomalous 
parameter used by the authors of [60]. Eqs. (7.36) and (7.38) are satisfied by our solutions 
(7.33) and (7.12) as long as one fixes the anomaly parameter A to zero. However Eq. (7.37) 
is satisfied with the gravitational anomaly switched on. So far these constraints are satisfied 
except Eq. (7.39), as ^5 is not vanishing in our model even though we fix the anomalous 
parameter to vanish. 

To check the value we get for ^5 (4.42), we may proceed by using the Kubo formula 
formalism. The Kubo formula for ,^5 will relate this coefficient to a two point function at 
second order in a frequency and momentum expansion. Actually it will appear in the same 
correlator as the chiral magnetic conductivity. To do so we can switch on a gauge field 
in the y direction Ay = Ay{t,z). In such a situation the Fourier transformed source J'^ 
reduces to 

j(5)" = iok.Ay , (7.40) 

so that using the constitutive relation we can read the two point function 

(J-JV) = -i^^k, + alEUjh , (7.41) 

where we have redefined Ie = ■^s/c for the reason we will explain below. We have checked 
that this Kubo formula leads to our result (4.42) for Ie by using the model of section 2 
considered in the probe limit. This is a non trivial check of the non vanishing of .^5. We 
leave a detailed analysis of this and other Kubo formulae for a forthcoming paper [69] . 

In order to understand the discrepancy between this result and the prediction done 
by the authors of [60], we can analyze the properties under time reversal of the source 
associated to .^5, which reads in the constitutive relations as 

jM = ^^ef^'P^.VpEx + ... . (7.42) 

This equation in the local rest frame u'^ = (1, 0, 0, 0) looks like 

J = |5Vx^ + ... = -|5^ + .... (7.43) 



21 



The electric field and the operator Vx are even under time reversal while the current is 
odd, in consequence the conductivity ^5 is T— odd. The fact that this transport coefficient 
is T— odd tells us that such a source might contribute to the entropy production. For 
this reason demanding a non contribution to the production of entropy might not be well 
motivated. The situation would be similar as demanding a vanishing contribution from 
the usual electric conductivity. One can see also the odd property of ^5 from Eq. (7.41), 
as {J^J^) is 7~— even and inverting the time is the same as changing a; — )• —u. 

We have noticed that the anomalous coefficients associated to sources constructed with 
the second derivative of the fields can be naturally factorized as 

Ai = -2riL , (7.44) 

A4 = -2r]lB , (7.45) 

Cs = (tIe . (7.46) 

These expressions make their dissipative nature clear, and they suggest the existence of 
anomalous relaxation lengths in analogy to the relaxation time Ttt- These new T— even 
quantities write 



/. 



|(^,,,,.(3.i_V-^)), (7.4, 



1 j 1^1-'' I 2rp2 ■ 



lE = -%{'^\og2-2\(l + 2\ogi))+0(p?). (7.49) 

A last interesting observation comes from the result on the dispersion relation of shear 
waves in [60], where they have found that 

oj7^-i^k^^iCk^ + ... , (7.50) 

4p 

with C = — Ai/(8p). It would be interesting to generalize the computation of [70] to the 
case including the mixed gauge-gravitational anomaly to verify whether the result for C is 

C = f L ■ (7.51) 

4p 

8. Discussion and conclusions 

We have studied the transport properties of a relativistic fluid affected by gauge and mixed 
gauge- gravitational anomalies. We have used a holographic bottom up model in 5 dim that 
implements both anomalies via gauge and mixed gauge-gravitational Chern-Simons terms. 
This model was used in [38] to compute the first order transport coefficients in holography 
from Kubo formulae. Within the fluid/gravity approach, in this work we have reproduced 
previous results at first order, and extended the computation up to second order in the 
hydrodynamical expansion. The computation has been performed in a Weyl covariant 
formalism, which allows us for a clear classification of terms contributing to second order. 
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We have found all the anomalous and non-anomalous transport coefficients of the 
model up to second order, except the ones associated to curvature sources. Most of the 
non-anomalous coefficients receive non trivial contributions coming from the anomaly sector 
through terms quadratic in the anomaly coefficients k. and A. There is a set of coefficients 
which are not affected by the presence of the anomalies. These are 



+A8S<''P^>/2 + A^E^^'E"^ + ... , (8.1) 

J^' = iicj^^^V^fi + i:iP^'^V''a^a. + C^^^^E^ + ^ju'^V^E'' + ... . (8.2) 



In particular Ai is usually redefined in term of the relaxation time r^r, in analogy with the 
the Israel and Stewart theory, Ai = —2r]TT^. 

On the other hand, we computed the anomalous coefficients which are (C, P) = (it, — ), 
i.e. these contributions are linear in k or A. These second order transport coefficients are 
dissipative unlike the ffist order ones, and in consequence some of them could contribute 
to entropy production. One example of that is the non vanishing value of .^5 (7.49). This 
coefficient was previously predicted to be zero by using non production of entropy argu- 
ments [60]. It would be interesting to compute within the present model the holographic 
entropy current and its divergence to study the contributions of such coefficients to the 
entropy production. 

We have defined the 7"— even quantities Ie, Ib and l^ in analogy with the definition 
of Ttt. a generalization of the Israel and Stewart theory to hydrodynamics with anomalies 
would give us a better physical intuition on these parameters. In particular the authors 
of [60] have noticed that the chiral coefficient C (7.50) in the dispersion relation of shear 
waves is related to Ai, and in consequence to l^, 



ojK. -ij-k"^ (l±Lk) . (8.3) 



The role played by the gravitational anomaly in the chiral vortical effect has shown up 
for the first time in the calculations of Kubo formulae involving the two point functions 
firstly derived in [71]. It is possible to write new Kubo formulae in terms of two point 
functions in order to compute the second order transport coefficients associated to second 
derivatives of the background fields. The rest of the transport coefficients at second order 
would demand the knowledge of three point functions. The advantage of Kubo formulae 
is that they can be used in either field theory or holographic computations. This would 
allow to compare the weak and strong coupling regimes of second order coefficients [69] . 
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A. First order solutions 

Here we show the exact form of the Fi[p\ functions defined in Eqs. (6.10) and (6.12) 

' (9Q«p^ + 2Af2p^(-4M + 3(1 + M)p) - 27Q^ (l + M - p^) + gQ^p (4m2 + 3^ _ 6Mp*)) p2 (2 + 5p^ + 2pf) 



Fi[p] 



8A/Vp2(-l + pi)(2 + pi)'(l + 2pi) 
IttQ (-1 + p2) (p2 ^ p2) (-1 + plf (1 + p2 + pI) (2 + Spj + 7p| + 5p6 + 2p|) 
8M2p6p2(-l+p2)(2 + p2)2(l + 2p2)=* 



3 (Q2 - Mp^ + p6) p2 (-2 - 3p2 + 3p6 + 2p|) ArcTan 

2Mp6(2 + p2)2(l + 2p2)3 
3Q 



/1+P2 



3(1 + A/)Q (Q2 _ Mp^ + pS) Log[l + p] 
2A./(-2 + Q2)2p6 



2Mp6 (2 + pIY (1 + 2p2)^ 2A/ (-2 + Q2)2 ^6 

(Q^ - A-/p^ + p6) (-1 + p2)p2 (2 + pI) Log[p - pa] 3Q (Q^ ^ A/p^ + p^) p2(l + P2) (2 + pj) Log[p + pa] 

4p6(l + P2)2(l + (-1 + P2)P2) (1 + 2p2)' 4p6(~l + P2)2 (l + P2 + pD (l + 2plf 

(q2 „ ^ip2 + p6) (-2 - p2 + A^p4) Log [1 + p2 + Pi] ^ ^ _, ^ 



3Q (Q2 - AJp2 + p6) (-2 -pI + Mpf) Log [l + p^ + pj 
4Afp6(2 + pi)2(l + 2pi)^ 



^ 2V3^p6(l + 3pi + 2p4)^ _ ^^^ / l2Afp| (1 + 3pi + 2p4)^ ^ 2Q2p2p2 (1 + 3^2 ^ 3^4 ^ 2p6)- 

'^''^ •3p6(l + 2pi)'(Q2 + Pl) V Qpio(l + 2p2)2(Q2+p6) 

+6p6 (Q2 + p6)2 (^ _^ 3^2 ^ 4p4 ^ 2p6) + AM^p'^pI (2 + 7pg + Qpf + 4p6 + 2p|) 

QplO(l + 2p2)^(Q2+p6) 



Op6(l + 2p2)2(Q2+p6) 

^ Qp"(l + 2p2)^(Q2+p6) 

~2p^pl (4 + Q^ (28 + Q^ (60 + (43 + 34Q^) pj (l + p|)))) \ 32V3A/^A (Q^ - Mp'^ + p^) Log[p] 

QpiMl + 2pi)'(Q2 + pl) )^ Q^' 

8V3QA (Q^ - Mp"^ + p6) (2 + 12p| + 27pf + 35pf + 27p| + 12p^° + 2pi^) Log[p - pa] 

P«p!(l + 2pi)' 
8V3QA (Q^ - A/p^ + p^) (2 + 12pi + 27p^ + 35p^ + 27p| + 12pi° + 2p^^) Log[p + p2] 

P«pni + 2pi)' 

8v^A (Q2 ^ A/p2 + p6) p2 (-1 - 3p| - 6p^ - 7p6 - 3p| + 2pf) Log [l + p^ + p|] 



p6(l+p2)3/2(l + 2pi) = 



(A.2) 



i^3[p] = V^ 



9Q5 + 2A'/Qp-* (-A/ + 2p (1 + 2p2)2) 
8Af2p6(i + 2p2)2 
(P2 + pi) (3A^Qp2p2 (2 + p2) - 27r (-1 + p2) (p2 _ p2) (1 + p2 + pi) (l + 4p2 + 6pt + 5p| + 2p|)) 



8A-/2p6(2 + pi)(l + 2p2) 
V3 (Q2 - A/p2 + p6) p2 (1 + plf ArcTan 



2A^p6 (2 + 5p2 + 2p4) 



p 

Vi+pI 



v^Q (Q^ - A/p^ + p6) Log[l + p] 
. V .-..■ ... ^ 2p6(-2-pi + 2p6+p|) 

V3Q (Q^ - Mp^ + p6) p2(4 + P2(-1 + 4p2))Log[p - P2] V3Q (Q^ - AJp^ + p6) pj (4 + p2 + 4pi) Log[p + P2] 

4p6 (1 + 2plf (1 + p3) 4p6 (1 + 2p2)^ (-1 + p3) 

V/3Q3 (Q2 _ Mp2 + p6) (_1 + p2)2 Log [1 + p2 + p2] 



4p6p2 (1 + 2p2)3 (2 + 3p2 + 3p4 + pI) 
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Fi [p\ = ^^ A 



Mp8(l + 2p2)2 \^ Mp8(l + 2p2) 

+2Q6p2 (_g7 _^ gp4) ^ 4 ^_5p2 _^ 3^6 + 6 (p2 + 2p^)^ (1 + 2pl + 2p4 + p6)^ 



,2n2 



Mp8 (1 + 2p2) 
24MA (02 _ Mp2 + p6) Log[p] 6Q2K/[p]Log[p - P2] 12A/[p] (l + p2) (1 + 5p2 _^ 9^4 _^ 5^6 + ^8) Log[p - pa] 



Q2p6 



(l + 2p2)^ 



p2(l + 2p2)^ 



6K/[p]Q2Log[p + p2] 12A/[p]Q2 (1 + 5p2 ^ 9^4 ^ 5^6 + ^8) L^gf^ _^ ^^j 

(l + 2p2)^^ p4(l + 2p2)« 

6g2K/[p]Log [1 + p2 + p2] 12A/[p]p2 (1 + 2p2 - p6 + p|) Log [1 + p2 + p2] 



(l + 2p2) 



2n3 



(l+p2)(l + 2p2)^ 



(A.4) 



^ / 3Qp2 (2 + 5p2 + 2p4) (-8M2p + 9(1 + M)p2 + 9(1 + M)pf, 

V 8M2p2p2(-l+p2)(2 + p2)2(l + 2p2)^ 



27r (3Q2 - 2Mp2) (-1 + p2)^ (2 + 5p2 + Tpf + 5p6 + 2p8) 
8M2p2p2(-l+p2)(2 + p2)2(l + 2p2)'^ 



\/3 (-3Q2 + 2Mp2) (-2 - 3p| + 3p^ + 2p|) ArcTan 



n+pi. 



\/3(l + M) (3Q2 - 2Mp2) Log[l + p] 
2M(-2 + Q2)2p2 



2Afp2yi + p2(2 + p2)2 (1 + 2p2)'' 
V^ (-3Q2 + 2Mp2) (-1 + p2)p2 (2 + p2) Log[p - P2] v^ (-3Q2 + 2Mp2) p2(l + p2) (2 + p2) Log[p + P2] 

4p2(l + p2)2(l + (-1 + p2)p2) (1 + 2p2)'' 4p2(-l + p2)2 (l + P2 + pf, (l + 2p2)=' 

V^ (3Q2 _ 2Mp2) (-2 - p2 + p4 + p6 + p8) Log [1 + p2 + p2] 



4Mp2(2 + p2)'=(l + 2p2) 



(A.5) 



^ 6Kpt(l + 3p2 + 2p4)^ ^ (*Q^ (120 + 77Q2 + 86Q4) p4 _ 4 ^_^pA + 9Mp4 (^ + 3^2 + 2p4)2) 



Mp2(l + p2)(p2+2p3)2 \^ Mp6(l + p2)(p2+2pi)=^ 

+8Q2p2 (^_28p2 + 3 (1 + 3p2 + 3p4 + 2p6)') \ 32M2A (3Q2 _ 2Mp2) Log[p] 



Mp6(l+p2)(p2+2p3)2 

8A (3Q2 _ 2Mp2) (2 + 12p2 + 27p^ + 35p| + 27p| + 12pi" + 2pi2) Log[p - P2] 

p2pt(l + 2p2)« 

8A (3Q2 - 2Mp2) (2 + 12p2 + 27p| + 35p| + 27p| + 12p^" + 2p^2) Log[p + ^2] 

p2pt(l + 2p2)3 
8A (-3Q2 + 2Afp2) (-1 - 3p2 - 6p4 - 7p6 - 3p8 + 2pi2) Log [l + p2 + p2] 
p2(l + p2)^(l + 2p2)3 



(A.6) 



Fs[p\ 



(1 + p2) (27p5 (1 + p2) (2 + pI) + 2-kQ (-3Q2 + 2Mp2) (l + 3p2 + 3p4 + 2p6)) 



8M2p2p2(2 + p2)(l + 2p2) 



,2\2 



n+pi 



Q3 (3Q2 - 2Mp2) ArcTan 



(3Q2_2Mp2)Log[l + p] 
' 2p2p3 (2 + 7p2 + 9p4 + 7p| + 2p|) "^ 2p2 (-2 - p2 + 2p6 + p|) 
(3Q2 _ 2Mp2) p2(4 + p2(-l + 4p2))Log[p - p2] (-3Q2 + 2Mp2) p2 (4 + p2 + 4p2) Log[p + pa] 



4p2(l + 2p2)^^(l+p3) 
(3Q2 _ 2Mp2) (-1 + p2)2 (1 + pD Log [1 + p2 + p2] 
4Mp2(2 + p2) (l + 2p2)^ 



4p2(l + 2p2)^(-l + p3) 



(A.7) 
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^'^''^ " p2(l + 2p2)2(Q2+p6) + Qp4(l + 2p2)2(Q2+p6) 

2y3QK (3Q2 _ 2Mp2) Log[p - pa] _ aVzQX (SQ^ _ 2Mp2) (l + 5pi + 9p^ + 5p^ + p|) Log[p - pa] 
^ pMl + 2pi? p2pi(l + 2pi)=^ 

2V3Qfv (3Q^ - 2Mp^) Log[p + pa] _ 4v^QA (3Q^ - 2Mp^) (l + 5p^ + 9p^ + 5p^ + p|) Log[p + pa] 
^ pMl + 2pi)' p2p!(l + 2pi)^ 

2 V3k (3Q2 _ 2Mp2) pa ,^1 + p^Log [l + p2 + p2] 4^;^ (3^2 _ 2Af p2) pa (l + 2p2 - pf + p|) Log [l + f?+ pi 

PMl + 2pi)' ^ p2(l+p2)3/2(l+2p2)3 

Sv^MA (-3Q2 + 2Mp2) Log[p] 
p2pi(l + p2)3/2 



{A.8) 



B. Second Order Sources 

In this appendix we show the second order sources in the vector and tensor sector. The 
are sphtted in terms of the anomalous and non anomalous ones. 

B.l Second Order Vector Sources 

In the vector sector the sources are organized as follow 

10 5 10 5 

^. = T.^i^^^i{'^ + Y.'^a^^^i^^ ' A, = j;rr)j-W+j:fr)jW, (B.l) 

a=l a=l a=l a=l 

where the tildes refer to the anomalous sector. The sources write: 
B.1.1 Non-anomalous vector sources 

3 (E) _ 2'r«T^(-l + p)^ (-p^l + p(2 + 3p)) + Q2(3 + 2p(3 + 2p(2 + p)))) 27,T^ p^ aQFi,[p] i7,T^ [-1 + p) {l + P + p'^) F[[p\ 

''+''^ ^ v^M(l + M)p8/[p]2 v^{l + M) v/3(1 + M)/[p1 

47rT2(-l + p)2 (Q2(g ^ p(;^2 + 7p(2 + p))) - p2(4 + p(8 + 3p{4 + p(3 + p{2 + p)))))) Fi [p] 

i '- , fB.2) 

v^(l + M)p7/[p]2 

^ 64\/37rQT2A2 (-SQ*" (242 + 131Q^) + MQ^ (2654 + 1603Q^) p^ - 4*1^ (l46 + \27qA p'^\ 

''+''2 " Af(l + M)pl7/[p] 

64x/37rQT2A2 (3 (16 - 498(5^ - 225Q'' + 16Q6) p8 + 220 (2 + 3Q^ + Q'') p**) 2v^7rT2p2Fi [p] A-nT^ p'^ F[[p] 
M(l + Af)p"/[p] (1 + M) \/3(l + M) 

TesySTrQ^T^^A (-3 + p2 + q2 (_a + p^)) fHV3TrT'^ \ (p^ + Q^ (-5 + p^)) agFrb] 16,rQT= ASqF^[p) 

M(l + M)pll/[p] (1 + Af)pf> (1 + Af) 

4S7rT2A (2p2 - 2p'3 + Q* (-5 + Sp^) - Q^ (s - Sp^ + p6)) F!,[p\ aTrQT^F^M ISttQT^ApSqF^'M 

Af(l + Af)p6/[p] (1 + Af)p (1 + Af) 

S^ttQT^A (-3Q^ (2 + Q^) + A/ (18 + 7Q^) p^ - 3 (l4 + 9Q^) p'^ + 8Afp**) F^[p] 16v^7rT^A (-5Q^ + 2A'/p2) aQF.J[p] 
Af(l + Af)pl"/[p] {l + A/)p4 

r+r<^) = p^F^Cp), (B.4) 

,E> 24192Q-»a2 - 76S0MQ2A2p2 + 12 96A (p^ + q2 ^5 + p2\\ ^^[^j a4A (-TQ^ + aAf p^) F^[p] 
r+r<^) = ^j-^ — + ^ ^ '-^ ^ -. '- + 32x/3QAF^[p] , (B.5) 

2 (E) _ v^'3(-l + P)^ (p^(l + P(2 + 3p))-Q2(3 + 2p{3 + 2p{2 + p)))) 2{-l + p) (l + p + p^) F^[p] 2{-l + p) (l + p + p^) F;^[p] 
''^''^ " Afp«/[p]2 J[p\ J[p\ 

, 2Q2(_i + p)2(5 + p(i2 + 7p(2 + p)))F3[p] 2p2{-l + p)2(4 + p(8 + 3p(4 + p(3 + p(2 + p)))))F3[p] 
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96V3QA (-999Q*^A + ISlBMQ^Ap^ - 978MQ'^Xp'^ + 2MQ'^(4k + 341X)p^ - 120M^Xp^^) ^ ^QF^[p] 
r.ri ' — ^^ — 2p' -Fo fpl + 

geySQA ('('l52A + 456Q^A + 152Q'^ \ - 3Q*'(4k + 121A)") p") 24A (qQ* - SMQ^p^ + M^ p* +3Q^p'' - M p^) F^[p] 
4\/3QA (gQ* - igMQ^p^ + GM^p* + 45Q^p'' - 22A'/p**) F^[p] 



3p"F3[p], (B.7) 



Af (-1 + p)p4(l + p) (Q2 - p2 (1 + p2)) 
ri.'^> = 0, {B.8) 

^ 3847rT^A (sQ* (l4 + Q^) A - 4MQ^ (22 + Q^) Ap^ - M^ (-18+ Q^) Ap^) 32x/37rQT^ A (p^ + Q^ (-5 + p^)) FgM 

+ * M(l + M)plS/[pl Af{l + M)p5 

3847rT2A (q^ ((e + 4Q2) k + 5 (i4 + Q^^ a) p^ + M (-18A + q2{-2k + A)) p**) 64\/37rT2A (p^ + Q^ (-5 + p^)) SQ-Fglp] 

M(l + M)pl5/[p] (1 + M)p5 

levrT^A (-Q2 (24 + 13Q^) + Af (e + TQ^) p^ - ((i + Q^) p'') F^lp] 

Af(l + A/)p6/[p] 
64,rQT2AaQFi[pl 2,rT2pF^[p] lex/SirT^A (-7Q2 + 2Af p^) 9qF^ [p] ISv^ttQ^t^a (-1 + 4p2 + 2p'' - Spi^ + 3p**) F^ [p] 

(l + Af) V^(l + Af) (1 + Af)p1 Af(l + A'/)plO/[p] 

16V3TrQT^\ (7p^ + p* - ISp*" + 3p** + Q* (4 - 3p^ + p^)) Fl,[p] 
Af(l + Af)plO/[p] ' 

3 (B) _ 48v/3QA (p2+g2 (^-5+ p2))F9[p] , 24a(p* -p** -Q%2 |-_g^p2^ (^_2^p4^ _l_p4|-j^_-^„p2^p4^^^^[pj 

8y3QA fl5Q* + 2Afp* f3Af - Tp*) + Q^ f-23Afp^ + 33p'^)) F^Ap] 

i i '- i '-^ — (B.9) 

AfplO/[p] 

o r,=.^ 8v^a2 ('-480Q3 +96Qp2 +96Q3p2) 96A f p^ + Q2p<l (-s + p^)) Fglp] , 8A f 21Q=^p5 - ep'' - GQ^p^) F/ [p] 
'■+'-$^> ^ 5 ^ + ^ ^^ ^-^ + 32x/3QAFi[p] + ^ 1 , 

p^ p^ pit 

(B.IO) 



3 (M) _ 'rT2(-l + p)(8p2 (l + p + p2) +Q'» (3-p(-3 + 4p(l + 2p)2))) ,rT2Q2(_i + p)(l2 + p(l2 + p(l3 + p + p2+9p3))) 
''+''^ ^ 2Af2(l + A<f)p8/[p] 2Af2(l + Af)p8/[p] 

^ 4^T^F^[p] 87rQT^(-l + p)(l + p + p^)Fi[p] 37rQT^p-V[p]3QF^ [p] 

\/3(l + Af)p2 (l + Af)p5/[pl 2A-/(1 + Af) 

3 (M) _ TT2(Q'5(-3 + 4p2)+8(p^+p6)+Q''(6 + p2(-l + p(-3-p + 8p3)))) 
''+''^ " 4Af2(l + A/JjpE' (Q2 - p2 (1 + p2)) 

7rT2Q2p2(14 + p(12 + p(2 + p(9 + p(-16 + 9p)))))) 47rQT2Fi[p] 327rQT2«;9QF7[p] irQT^ (3Q2 - TAfp^ + 3p'') F^ [p] 
4A/2(l + Af)p5 (q2 _ p2 (1 _|_ p2)) (1 + A<f)p2 (1 + A'/)p3 2A<f(l + Af)p3 

24v/3,rQT2A(-4-3Q2+2Afp2)F^[p] lev^vrT^A (-5Q=^ + 2Afp2) SqF^[p] 37rQT2p3/[p]aQF^[p] 3v/37rQ2T2/[p]F^ [p] 

A/(l + Af)p<l (1 + A/)p«l 4A/(1 + M) 2M(1 + M)p 

8irT^ (-6Q^ (4 + ■iQ^) A + 12AfQ^Ap2 - f 2 + Q^) Kp'^) F'[p] 3S4TrQ^T^\^ f-4 - SQ^ + 2Mp^) f-63Q^ + 20A/p2) 

+ i i '- i '- '- + i '-^ '- ,(B.12) 

Af (1 + A'f)p8 Af (1 + Af)pl= 

(M) 3V3Q 

r+i~\ = T. (B.13) 

^ -■ 2Afp2 

\/3Q (p^^ - 3840A2p* (-1 + p'^^^ 48\/3QA (-5Q^ + 2Mp'^\ F^[l\ 16\/3Q (-ISQ^A + SA'fAp^ - Kp") F7[p] 
pl5 p9 p9 

24A (sQ^ - 2A'fp^) F^[l] \/3Q (384Q^A^p^ (-46 + 20p^ + 25p* + p^') + 192Q''a^ (i()5 - 92p^ + 20p^ - 55p'' + 22p**)) 

(-l+p)(8p'' (l+p + p^) +Q''(-27 + p(-27 + 4p(l + p + p2) (-9 + 2p2)))) 2 (-1 + p^) F^[p] 

4 Af 2 p8/[p] -^^^f"' Ji 

(-1 + P) (Q^P^(-9 + p{-9 + p(7 + p{43 + 16p))))) 4\/3Q{-l + p) (l + p + p^) F3M 

i '- i '- , (B.15) 

4Af2p8/[p] pS/[pl 
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(M) 
'•+'•4 



(B.14) 



3Q2 (512MA2 (-3Q2 + 2Mp^) (-63Q2 + 20Mp2) + gQ^pl" - 6Mpl2 _ 9^13-) 3v^Qi?^[i] 9Q'^ Fs[l] 

8 M2pl6 8Afp2 4Mp2 

2v^QF3[p1 1 12v/3QA(-3Q2 + 2Afp2)F^[p] v^Q (3Q2 - 7Mp2 + 3p6) F^[p] 

7^ 2^^'"' + 177 + S?;^ 

4 (-3Q2 + 2Mp2) [(>Q^\ + «p6) F^[p] (-4p** + 9Q'' (-1 + p^) + Q^ p^ (9 - 9p* - 4p«)) Fl,[p] 



Mp^ AMp'' 



{B.16) 



^, , ^ _ 37351^ _ 9|^ _ ^ _ ^ 
4p2 2Mp2 

^ 3v^7rQT2 (vrTplI +512r+A2 (35Q* + 15Q'' +2M2p=^ (-2 + Afp2) +q2 (20 - 13M2p2))) ttT^ (9Q2 + 2Mp2) Fi [p] 

+ * ^ 2r+A'/2{l + M)pl3 \/3bO«'M(l + A<f)p2 

le^rQ^T^^FoM ■A27,QT'^naQFc,[p] 2^T^ pFilp] ■■i2V37rQT^ X (^-3 + p^ + Q" [-1 + p^)) Fi[p] 
M(l + M)p3 (1 + M)p^ \/3(l + M) M(l + Af)p«l 

WVSTyT^X (-5Q^ + 2MpA aQF^[p] SttQT^ (q^ - Mp'^ + p*^) F^lp] 

+ ^ \ ^ , (B.IB) 

(1 + M)p'i 2M(1 + M)p'^ 

3 (M) _ 3v/3Q(i024Ma2 (l5Q* -13MQ2p2^2Af2p4) +9Q2pl0 _2Afpl2) ^ (gQ^ + 2Afp2) FgM 24,Q'^KFg[p] 

'■+''!' = SIJV^^ ''''^'''' 2]V?;^5 + ^ 

48V3QA (p2 + q2 (-2 + p2)) F^[p] 3v^Qp/[plF^ [p] 4 (-SQ^ + 2Afp2) (eQ^A + ,^p«) F^ [p] 

H — -|- , (B.19) 

A'/p<l 4Af A'/p** 

2 iM) _ 1 16v/3Q^F9[p] 24A(2p^+Q^(-5 + 2p^))Fi[p] 
'■+'"10 - ^T; 73 Ti ' (B.2U) 



B.1.2 Anomalous vector sources 



-(E) 
'•+'•1 = 



16x/3QA(2p2(-3 + p(-3 + p(-3 + 2p))) + Q2(17 + p{17 + p{ll - 10p{l + p))))) 16v^QA (-63(3^ + 20A<fp2) F5 [p] 

pH-l + P) {-Q^ + P^ + P'') 7 

2p= (-Q^{6 + p(12 + 7p(2 + p))) + p2(4 + p{8 + 3p(4 + p(3 + p(2 + p)))))) F2 [p] 2(-l + p) (l + p + p^) Ff,[p] 

(l + p)M-0"+P"+P*)' ?W ' 

(B.21) 

2 _(g) _ 48A(p2+p3+p4 _p5 + Q2(_4 + p(_4 + p(_3 + 2p(l + p))))) 96A (p^ + Q^ (_5 + p2)) F5 [p] 

'■+''2 " p3(l + p)(-Q=^+p2+p4) ^5 

2p= (-Q^{6 + p(12 + 7p(2 + p))) + p2(4 + p(8 + 3p(4 + p(3 + p(2 + p)))))) F4[pl 2(-l + p) (l + p + p^) F^M 

(B.22) 

^ ,„, -120Q^A+ v^QfgM +48Ap'' +48Q^Ap^ 
'■+4 = ^ e 3p^F4[pl - 2p3F^[p] - pF;Ip] , (B.23) 

^ ,^, Sv^ttQT^A (2p2 (-5 + p2) +Q2 |^_21 + 7p2 +2p'')) ISV^ttT^A (9Q* - BAfQ^p^ + A/2p<l + 3Q2p6 _ jy-p8~J j^^ [pj 

''+'''' ^ Af(l + A/)p8/[p] Af(l + A'/)p6/[p] 

(i4s/37TT^\ (p2 + q2 (^_5 + p2^-) SQFgM 32v/3,rQT2A (p^ + Q^ (-5 + p^)) Fglp] fATrQT^ XaqF^^lp] 

(1 + Afjp-'i Af (1 + Af)p= (1 + M) 

leirT^A (q* (13 - 7p2) + 6p2 (-1 + p*) + Q^ (24 - 13p2 + p")) F^ [p] ISv/SttT^A (-7Q2 + 2Afp2) SQ-Fglp] 

Af(l + A/)p6/[p] (1 + Af)p4 

SttQT^A (gg* - igA/g^p^ + eA^^p* +45Q^p'' - 22A'/p**) F|![p1 ISv^ttQ^T^A (-1 + 4p2 + 2p* - Sp'^ + 3p^) Fl^[p] 

Af(l + Af)plO/[p] Af(l + A/)plO/[p] 

16\/37rQT^A (7p^ + p* - lOp" + 3p** + Q* (4 - 3p^ + p*)) F|^[p] 



Af(l + Af)plO/[p] 
96A(p2 + Q2 (-5 + p2))F8[p] 24A(2p2+Q2 (_7 + 2p2) ) F^ [p] 



(B.24) 



i-^?^ ' = 32V3QXF^[p] + 

(B.25) 



28 



_^^j _ 12A(-63Q*p + 8M2p2+2MQ2 (^-8 + 15p3)) 4\/3Q(-l + p) (l + p + p^) F2 [p] 3^Qp^F!^[p\ 2 (l - 2p3) F^[p] 
"'' ^ ^ Mp9 P^/[p1 2M p2 

48 (i5Q*A + 4A'/2Ap*' +Q^ (-16MAp^ + Kp'')) -F5 [p] (sQ^ + 2Mp2) Fjlp] 
^ ^;il Mp2 ' 

2 ,(M) _ 4v/3Q(-l + p)(l + p + p2)F4[p] l6v^Q«F5[p] 3x/3Qp3F^[p] 2 (l - 2p3) F^ [p] 
'■+''2 " ^?7W ^ ^^ 2M + ^5 

36v^QA (-SQ^ + 2MpA - hQ'^ p'^ + 2Mp'^\ Fg[p] 
Mp<^ 

2 -(M) 9Q^-F2[1] 9^/^Q^F^[l\ 3v/3Fi[l]Q 9Q^F9[1] / 9<3M „ r , 'iVSQFi[p] 9^Q^ F-j[p] 

-"-"—— — ^ ^ -^'1+ ^2 [Pi — 



4M2p2 2M2p2 2Afp2 Mp2 Y 2Mp'^ J p2 4M2p2 

gySQA (iOMQ2 (q2 - 2p'5) + p* (4 + 12Q2 + 33Q^ + 4Q<' - 2 (4 - 12Q^ + 5Q^) p")) 

M2pl2 
9v^QAm2 (-14Q2p2_|_4pS^ QQ^FgM P (8Af2 + 9Q2p2^ ^^[^j V^Q (-3Q^ + TMp^ + 3p'') F^[p] 

Af2pl2 Afp2 gM2 4Mp3 

3-/3Q (q2 - Mp2 + p6) F^[p] 9Q2/[p]j7^[p] 

4Mp^ 4Mp 

^18^^Q^t3x(^-21Q'^+10Mp^) tt-^T^ (9Q^+2Mp^)Frlp] IQ^Q^T^ ^Fglp] 327rQT^ ^dQFslp] 
r+M3{l + M)p8 2r_|_Af3(l + M)p2 M(l + M)p3 (1 + M)p3 

24,rQT2A (-3Q2 + 2Mp2) F[[p] 3V3tt^ QT^' p^ F^lp] ■i2s/3nQT^ \ (-3 + p^ + q2 (-1 + p2)) f;,[p] 

M(l + M)p''- 4r+Af3(l + Af) (1 + Af)p* 

leVSnT^X (-5Q2 + 2Afp2) dqFl^lp] SttT^ (-3Q2 + 2Afp2) (sQ^A + np^) F^[p] 2^^T^pF^lp] 



(B.26) 



(B.27) 



(B.28) 



{B.29) 



(1 + Af )p* V3M(1 + M)p^ r+M2(l + Af) 

87rT^ (6Q^ (4 + 3Q^) A - 12MQ^Ap2 + (2 + Q^) Kpl^) F,^ [p] 
M(l + M)p^ 

9Q2f^[1] _ dVaQ'^Frll] 3-/3QFill] 9Q^Fgll] _ 9V3Q (4 - 12Q^ + 5Q*) A 
8Af2p2 4Af2p2 4Afp2 2Afp2 A<f2p2 

V3QFjIp] 16s/3QkFs[p] 24A(-5Q='+2Afp2)F;^[p] V3QpF^[p] 
— -^ +Folp] -3 +2pFg|p]. {B.30) 



B.2 Second Order Tensorial Sources 



In the tensor sector the sources are also sphtted in terms of the anomalous and non anoma- 
lous, the tildes refer to the anomalous sector 
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E^'^^i^^ + E^'^^i"^- (B.31) 



a=l a=l 
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B.2.1 Non-anomalous tensorial sources 



Pl = 2p-3p^F5[p] -2pV^[p], (B.32) 

(Q2-Afp2+pi3)F^[p]2 

P2 = 4p - 6p^F5[p] + ^ '- , (B.33) 

P 

76SAQ* fp^ f-SMA - 3{k - 29A)p'') +2Q2a ('9 + 4p'')) 32\/3QA ("q^ - 2p<') F;^ [p] 

^^ = pii(Q^-Mp2 + pO) ^ + ^-'^-^["l + ^-'^^W + Q2l,,,.Je 

2Q2 ('9216A2(p2 - M^) + 768M(k - 10A)Ap* + p**) + p^ ('3072A^ - M p^ + p^^) 32A (-3(3^ + 2Mp2) (q^ - 2p'') F^[p] 
~^ p7 (Q2 - Mp2 + p6) + pi (Q2 _ Mp2 + p6) ' 

{B.34) 
2/ 4/644352M 396288 \ g/ 549504 3072 \ 3072 \ 2Q^ M 22/ 130560M^ 9216 107520M 

, ^ , , , 2\ 16A(27Q2+4Mp2)/[p]F^[p] 
+ /[p] (i92V^QAF^[p1 - p-V^[p]2j + i '- P^F^lpf , (B.35) 

.^n = -^^|^(pV,[p])'. (B.36) 

/ ^T \ S,rQT3(p3Fi[p])' 3Qp2Fi[p] ., w , , , v^ (2 + Q^) F^ [p] 

y r_|_(l + M)y 3r;?(l + M)^ M ^ ^ Mp 

r+P7 = -2(p^F3[p])' , (B.38) 

3 6wQT^p^Fi[p] 4v^7rQT2p2F3[p] 4V37:T^ p^BqFsIp] 2nQT^ p^ f;^[p] 

r I Ps — — ^ ~ ^ ~ — 7= (B.39) 

+ r^M(l + M) M(l + M) (1 + M) V3M(1 + M) 

4,rT2p5F^[p]F^[p] ^ 4wT2p39QF^[p] ^ 4^T2p--i J[p]F^ [p]F^ [p] ^ 2^ (2 + Q^) T^F^ [p] 
V^(l + M) V^{1 + Af) v^(l + M) Af(l + M)p 

2 6v/3Qp^F3[pl , V^Qp^F^M 5^,, ,2 , (3Q^-2aV)f^[p] 32 ,„ ^„, 

r+Pg = + P F3[pl + — p flp]F,lp] , (B.40) 

A^ ('-57600Q* + 59904MQ^p^ - 8448M^p'' - 34560Q^p'' + 5376Afp**) 1 768Q2reA 

r.Pio = ^ TT, + z +48Ap2/[plF^[p] - p'^F;^[pf 

p^^ p p' 

112\/3QA/[plFo[pl 3 , , 

+ .,' ^ '^ ^ - p'/[plF^[p]=^ , (B.41) 

P'' 

256\/3QA ('3OQIA - -iHMQ^Xp^ + tiM^ \p* - 3Q^(k - lOAjp" + 2M{k - 3A)p**") 32v^QA (q^ - 2p'') Fl[p] 

""^^^ " ?f=7w ^«7w 

32A (-3Q^ + 2MpA (q^ - 2p<^) F(,lp] 

i i-^ i , (B.42) 

bOplO/[p] 

2v^QA^ f 112704Q'' - 127296AfQ^p^ + 24320A<f^p* + 71232Q^p'' - 16128Afp^") 2\/3QA f576Q^«;p'' + 128Afftp*) 

'•-1-P12 = ^ ^ 

+ ^^ pl5 pl5 

2s/3Q 2 , 2(96\/3Qa(q2 - Afp2+p6) -pllF;^[pl)F^[p] 112v^QA/[p]F' [p] 
+ ^=^ +48Ap"/[plF;5[pl + ^ ^ ^ ^ + ^ '"' ^'" 

pJ pD p^ 

16X (27Q^ +4Mp^) f[p]Fl,[p] 
+ ^ 7 2p^/[plF;[p]F^[p] , (B.43) 
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B.2.2 Anomalous tensorial sources 

Pi = g 2(p-3F2[p]) , (B.44) 

2_ _ 32,rT2A(-,rT(51Q2 +4Mp2) +i2r+MQ(Fi[l] + 2V3F6[l]Q)p3/[p]) 2^,rQT2p2F2 [p] 4^,rT2p2sQF2 [p] 
""+ ^ r_|_M{l + M)p» M(l + A/) (1 + M) 

768^7rQ2T2A/[p]F6[p] ^^T^ (2F{ [1]M + 4V3F6[1)MQ + 5Qp^) F^ [p] St^QT^ FeMF^lp] 

(1 + M)p^ V3M (1 + M) (1 + A/) 

4,rT2p3eQF^[p] 32,rT2A(27Q2+4Mp2)/[p]F^[p] iTr^T^^ F^lp] 4,rT%3 / [p] F^ [p] F^ [p] 2,rQT2p*F;5' [p] 

V3(l + M) •/3(1 + M)p'i r^M (1 + M) p VSil + M) V3M {1 + M) ' 

2. 16-^Q'^T^x(-1+p^ (3-5p + 7p'> -ip^^^ 647rT2A(2 + p* (-3 + p^ (-5 + 5P + P*))) 

""+ ^ ^ M(l + M)pl7/[p]2 M{1 + Af)pl3/[p]2 

le-irQ^T^A (-6 + p2 ("13 + p (-6 + p (s + p - 33p2 - 8p^ + 14p^ + 39p-°^ - 22p^)))) 
^ Af(l + Af)pl7/[p]2 

32TrQ2T2^|-5^p2 (^8 + p (-3 + p (-23 + p (-2 + p + 29p2 - 21p^ + dp'' +P**))))) 
^ M{1 + Af)pl5/[p]2 

128^7rT2AQ2 (q2 _ jp") Fglp] 327rQT2A (-5Q2 + 3A/p2 + p^) SqF^[p] 

(1 + Af)pll/[p] (l + A/)p6 

87rT2A (4p'' + 3Q* (-5 + 4p2) + 2Q2 (_5 + 6p2 _,_ ^6^^ ^^ j^j (UttT^X (3Q2 - 2Afp2) (q2 _ jpB^ ^^ [p] 

M(l + Af)p6 \/3(l + A'/)plO/[p] 

64F;[l],rT2AQ (q2 _ 2p6) 128y3F6[l],rT2AQ2 (q2 _ 3^6-) 32^Qt2 Ap/[p]aQF^'[p] 
(l + Af)pll/[p] {l + Af)pll/[p] ^ (1 + A/) ' 

r+P4 = -2(p^F4[p])' , (B.47) 

3- _ 16y37rQT2A(-18- 11Q2 +4Mp2) 4^7rQT2p2F4[p] 4^7rT2p2aQF4[p] 32^7rT2 Ap2/[p]F^ [p] 27rgT2p3F^[p] 
""+ ^ ^ A/(l + Af )p6 M (1 + A/) (1 + M) (1 + M) v^M (1 + A/) 

4,rT2p5F{[p]Fi[p] 4,rT2 p^agF^M 224,rQT2A/[p]F^[pl 2,r (2 + q2Jt2f^[p] 4,rT2p3/[p]F^[p]F^[p] 

— -|- -|- — , (t>.4ol 

^(l + Af) V3(l + M) (1 + Af)p2 Af(l + Af)p \/3(l + A//) 

56y3QA(-3Q2+2Afp2) 6y3gp2j,^;p, ^ ^^ ,2p,,,U',l 

r , Pe = r— s ^ H — ^ 48Ap flplF^lp] + p —— 2p F^[p] F^\p] 

M p^ M \ M J 

112V3QA/[p]F^[p] (3Q^ - 2A//p2) F;,[p] 



-V/W-FgW-FgW, {B.49) 

p- ivi p 

?,2VSQ\{q'^ -2p'^'] f!^[p] 12\(ip^ +Q'^ 
r+Pi = 



p2 Mp 

32V^QA (q2 - 2p'') Fl^[p] 12A (4p'* + Q'' (-5 + 2p2) + 2Q2 (p2 _ p6 ^ jpS^J^ ^^ [pj 



p6/[p] A/p6 

32a(-3Q2 +2Afp2) (q2 _2p6) F^[p\ 12A (-4p8 +(3* (_5 + 4p2) + 4Q2p2 (1 + p* -p")) F^' [p] 

4a(-3Q2 +2Afp2) (-18Q2p2 - 4Af p* + (q2 _ ^,jp2 + p'^'Y Fi^'> [p]\ 
A/plO/[p] ' 



(B.50) 



3^Qp2F2[pl r r- r , s , (-8A (27Q2 + 4Afp2) +p7F^[p]) (-3Q2+2A/p2+2Afp''/[p]F^[p]) 
r+Fs = „ - 2 (-96n/3QA/[p] +p-F^[p]) F^lp] - ^ ^ 1 L^ ^^ 



{B.51) 



C. Transport coefRcients at second order 

In this appendix we will write the expressions for transport coefficients up to second order. 

C.l Vector sector 

The solutions for the non anomalous coefficients ^i , . . . ,^io as written in (7.1)-(7.10), are 
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given in terms of functions ^ito^K'^^KX,x^) whose expressions are 



S.2,K\{p2) 



2nT-^ 



5GM^M + 1)Q4(1 + 2p1)^rl 



2430 - 14121M + 32625M^ - 36279M^ + ITISIM* 



)1 DA/f^n -I- 2o^)'^ 
+ n2 (^ + 6M - 7M^ + 2M^) log[l - pi] 

_^ 10M_£2(24_P2) /g^^ _ ^2^ ^ 3^^g _|_ ^3^2^ _ g^2^^3 _ ^^2^ ^ 8M^(17 - Upl) 



C2,A2(P2) 



Q2(i + 2pi) 



-32M*(4 - Pa) + 48M'' log 



nr 



2 + pi 



1-Pi 



(C.l) 



210GM3(M + l)p|(l + 2pi)4r^ \ M'^ + Q^ + pi(2M - 1) 



1837080 + 15089571M 



-54047817M^ + 109739475M''* - 13661086571^" + 102222345M^ - 35693475M® 



-7547847M' + 12206741Af" - 3911944]^^ + 427856M"' + 



o\ , 1680(1 -pi) 



Pi(l + 2pi) 



81(1 + 2p^ 



-M(351 + 783p2) + M^(540 + 1674pi) + M^(225 - 1719pi) - M'*(1497 - 672p2) 



+M'^(1332 - 12pi) - M'^(340 - 32pi) + 32M^) log[l - pi] - 



1680MVi(2 + pi) 
(l + pi)5(l + 2pi) 



X (81(1 + 2pi) - 27M(16 + 29pi) + 54M^(21 + 31pi) - 9M^(216 + 191pi 



C4.fcA(p2) 



+3M"(723 + 224p2) - 12M^(112 + pi) + 4M^{93 + 8pl) - 32M^) log [2 + pi] 
2^3 



(C.2) 



57rGM2Q3(l + 2p 



2V3 



(1 + 2pi)(405 - 1944M + 3708M^ - 3612M'^ + 190371^" - 507M'^ 



+57^4") + 



Q. 5M^{M'^-3){l + 2plf 



Q^pI 



<-3„6/r, , 2\3 



^ 5Mp2g+P2) ( 3^2 _ 3^(^ ^ 2^2) ^ 2^2 ) j^g 



2 + p^ - 5M + M'(3 - 2p^) log[l - pi] 



2 + pi 



C4,A2(P2) = 



1-pi. , 

1224720 - 8421651M + 24887169M^ - 41179203Af^ 



(C.3) 



280^/37rG(57M2(l + 2pi) 

+41512749M'' - 25998933M^ + 10023123M® - 236U97M^ + 330895M* - 21092M^) 

_ 3360M"(l + 2pi)" ^Q - 60M + W5M^ - 70M^ + 20M" - 2M^) log[l - pi] 



3360M''p|(2 + pi)^(Mpi - 3Q^ 



Q2(l + 2pi) 



-(-3pi + 3M(1 + 2pi) - 2M^) log 



2 + pi ^ 

1-pi 



(C.4) 



C5,o(P2 



3Q2 



2567rGM3(l + 2pi 

2MVi(2 + pi)^ 
3g2(Q2+pf)(l + 2pi) 



^ , (5M - 3)(11M - 15) + 8^'a + 2P^)' iog[i _ pi] 



(9pi + M(l - 16pi)) lo; 





3Q2 




r2 + pii 


s 


[i-piJ 



(C.5) 



C6,o(P2) = 



5127rG(l + Q2)3(l + 2pi)2 



'243 - 513M + 189M^ - 39M^ + 152M'') 



+8M^(l + 2pi)Mog[l-pi] 



2M^ 



(l + 2pi: 



■ (3pi + M(4 + pi)) (15(1 + pi) - 16M) log 



2 + pi 



1-pi 



(C.6) 
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£,6,K.\{p2) = 



107rGM3(52(l + 2p2)2 



+ 



2430 - 9477M + 13824M^ - gSTeAf" + 3444Af* - SOlAf' - 14M'' 



.2 ,.3^,_^, „2, , 10MVI(2 + pif 



Q^ 



(18 - 27M + 12M^ - M^) log[l - pa] + 



0^(1 + 2pi) 



2 + pi 

1-pi 



C6,a2(P2) 



612360 - 3485349M + 7477191M^ - 6516837Af'' - 444309M^ 



(3p2 - 3M(1 + 2p2) + 2M^) (3 + 6p2 - M(4 + Spa)) lo; 

1 

"2807rGM3Q6(l + 2p2)2 

+5355933M^ - 4409403Af'* + 1747097M'^ - 359455M* + 29492M'' ) + \, ^^'^ x 

1680MV^"(2 + pi)^ 
Q2(l + 2pi) 



(C.7) 



C7,o(P2 



(81 - 315M + 453M^ - 303M^ + lOlM^ - 16M^ + M®) log[l - pi] 

[ipl - 3M(1 + 2pi) + 2M') (^9(1 + 2pi) - 12M(1 + 2pi) + 3M'(1 + pi) - M^ ) lo; 

7rGMa(l + 2 ^X^ ( ^^^ + ^^^^' + ^^^^^ + ^^^^' + 128Q«)(1 + 2p^) 



2 + pi 

1-pi 



256 



+M^ (90 - 168M + 82M^ + 4M(1 + 2pi)^) log 



2 + pi 

1-pi 



(C.8) 



(C.9) 



C8,«;2(P2) 



^/37^Q^^ 



(,S,K.\{p2) ~ 



2GM3(l + M)(l + 2pi)-V3 
+2M(3 - M)(3 -M + 4M^) lo; 



3 , (18 - 69M + 63Af^ - 28Af^)(l + 2pi) 



2 + pi 



1-pi 



(C.IO) 



^».,.-, ,.N^r,, „ .TTT-^I Q^(216-1233M + 2691M^-2823M=' + 1379M'' 
GMi{l + M)Q->(1 + 2p^)5r| ' ^ 

-202M'^ - 24M®] (1 + 2pi) + 2(3 - 4M)^M^(1 + 2M - M^)(l + 2pi) log[l - pi] 

'M(216 - 1062M + 2094Af^ - 1950M^ + 634Af* + 264Af^ - 244Af'^ + 48M'^ 

2 + pi "' 

1-pi 



-18M^pi + 12M^pi + 82M^pi - 112MVi + 32M''pi ) log 



(C.ll) 



C8,a2(P2) 



ttT" 



35^/3GM3(l + M)Q7(i + 2pi)4r 



102060 + 846855M - 2955555M^ + 5599602M^ 



-5989080M* + 3085503M'' + 60033M® - 797012M^ + 279586M* - 32072M^ 

+ ^o^^^^ „, (3 - 4M)^(18 - 93M + 207M^ - 191M^ + 54M'^ - 6M^ + M^) log[l - pi] 
g^(l + 2p^) 

420Af / •) o o ■) o 

+ ^,,^ . „ „, ( 243 + 162M(-11 + pi) - 27M^(-207 + 47pi) - M'X9243 - 4383pi) 



Q2(l + 2pi) 

15M'*(-497 + 545pi) - M^(372 - 8382pi) - 2M''(2243 + 2203pi) + 9M^(404 + 113pi 



-M*(1219 + 120pi) + M'^(191 + 16pi) - 12M^'' ) log 



2 + pi 

1-pi 



(C.12) 
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i9,K^{p2) = 



-%/3Q 



4.GM3(1 + 2plf ' 9(3 - ^^^ + ^^') + 3(3 - mM' Iog[l - p^] 
"2(lTV) ' 54 - 9A/(10 + 2p^) + 2AM\2 + pi) + AM') ] log 



2 + pi "' 
1-pi 



(C.13) 



6,ka(P2) 



2%/3 



M3(-3 + 4M)Q 



[648 - 2079M + 2475M^ - 1365M^ + SSTM"* - 52M^) 



+ ^;,/, •„ ,,,■' ('2(1 - pi)^(l + Pa) ( - 36 + 129M - 129M^ + 28M^ 



M(-3 + 4M) 

+P2(36 - 39M - 12A/^ + 8Af^)) + 2M(324 - 1269M + 1935M^ - 1443M^ 
+537M'' - 92M^ + SAf** + piA/(-3 + 4Af)(-18 + 9Af + Af^))) log[l - pi] 

+ ^^, ~" ^,,, f - 972 + 5103A'/ - lOSSlA/^ + 12069Af ^ - 7383Af'' + 2424Af'' 



6,A2(P2 



2Af 

-392A/'' + 32Af ' + paAf (3 - 4A/")^(-18 + 9Ar + Af) ) log 
\/3pi 



'>) 



2 + pi 



(C.14) 



1407rGA4'3Q7(i^_2pi)3 
2Af^ + 51 
420A/2(l + 2pi) 



(-1 + 2A/ + pi) 102060 - 496935A// + 963855Af^ 



954012Af^ + 518949Af'' - 159036A//^ + 27094Af ^ - 2411Af '^ + 1&M^\ 



Pi 



243 - 1269Af + 2673A/^ - 2919A/^ + 1779A/'' - 610Af'' 



+11QM^ - 8AfM log[l - pi] + 

420A/" / T -) O Q T 

+ — 2— (243 - 243Af (6 - pi) + 27Af^(134 - 47pi) - 9Af^(529 - 297pi) 
P2 ^ 

+Af'*(3513 - 2919pi) - 3Af^(456 - 593pi) + Af'^(181 - 610pi) + 10A//"^(5 + llpi 



-4Af**(5 + 2pi) + 2A.fniog 



2 + pi 



1-pi 



CiO,k2(P2) = 



9pi 

47rGA/2(l + 2pi)3 
-|Af'(3 + 6pi + Af(5 + pi))lo: 



(3 - AM + 3Af^)(-l + 2A-f + pi) - 
"2 + pi 



2, 2Af2(l + 2pi)3 

i 



3pi 



1-Pi 



CiO,ka(P2) = 



27rGAf2(l + 2pi)2 



108 + 333Af - 372M ^ + 163Af^ - 16Af'*) 



(C.15) 
log[l - pi] 

(C.16) 
(C.17) 



+6(3 - 4Af )Af " log[l - pi] + 



2, 6Af2(6-5A/-A/2 + (3-4Af)pi) 



l + 2pi 



log 



2 + pi "' 
1-pi 



ClO,A2(p2) 



AQ^ 



35A'/2(i + 2p2)2 



51030 - 308205Ar + 782505Af^ - 1075848Af^ + 866631A/'' 



\ 42nA/f"^Cl — r?\'^('\ + o2~)4 

-421824Af ^ + 126626Af '' - 22369Af^ + 1874A/'** ) + \, ^^^ ,^ 

/ pi(l + 2pi) 

X (9 - 3Af (7 + pi) + 6A/'(2 + pi) + M^) log[l - pi] + ^^^^^g,^ (27 + 54pi (C.18) 



and for the anomafous coefficients ^1 , . . . ,^5, (7.11)-(7.15), one has functfons S,i.{K,x) that 



34 



write 



Ci,a(P2 



2ttT^ 



GM2Q2(l + 2p2)2)r2 



I ^9 _ 24M + 14M' + SM^'j + -—{1 + 2pl) log[l - pi] 



M^pl 



C2,a(P2 



(l + p2)(i + 2p2)2 
2GNPQ{l + 2plYrl 



(3pi - 3M(1 + 2pi) + 2M2) log 



2 + p^ 



1-p^ 



(C.19) 



2M' 



( - 6 + 7M + M^) + -—(1 + 2p^)2 log[l - pi] 



2Mpt 



^iAp-i) = 






(Af2 - 30'p^) log 



2 + pi 



SGM*{M + 1)(1 + 2pl)-irl 
2M'^{M -3) 



(5M - 3) (15 - UM + 4M^) 



(C.20) 



3Q2(i + 2p; 



(3-16M + 12M^)log 



2 + pi"! 



1-p^ 



(C.21) 



C4,a(P2 



4GM4(M + l)Q2(i + 2pl)'^'rl 

^6\ 27rM3(M-3)(2M + l)(l + 2pi)* 



+ 



1096M° + lOOM") 
27r(M - 3)M^pi 



1215 - 6075M + 11880M^ - 11367Ar + 5355M* 
log[l - pi] 



^5AP2) = 
isAPi) = 



4GM3(l + 2p2)2r2 



(9(1 + pi) + 3M(1 + 6P2) - 14Ar (1 + 2pi) + 4M^) log 



2 + pi 

1-pi 



2M2 



302(1 + 2p2 



log 



2 + pi 



I-P2 



^ttT^ 



2GM3(5(1 + 2p2)2r:^_ 
2M=*(4M-3)(0^-pi) 



^^^-^ I ,^27 - 72M + 54M^ - 7M^) 



+ 



QVi 



, ^^ 2, 2M3pi(Mpi-3Q^) , 
^°g[^'^^J' Q^(l + 2pi) ^°g 



2 + pi 



1-pi 



(C.22) 



(C.23) 



(C.24) 



C.2 Tensor sector 



In this sector the non anomalous coefficients Ai , . . . , A12 written in (7.16)-(7.27) and the 
anomalous ones Ai , . . . , Ag (7.28)-(7.35), are given in terms of functions ^i^[o^K2f^x,\'2) ^^^ 



A, 



i,{K,\) 



respectively. The expressions for these functions are in general very complicated, 



and we present here the result as an expansion at footnotesize p2 up to order 0{p2), which 
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is equivalent to order 0{ix^). For the non anomalous coefficients we get 
A4,;,2 (p2) = j^ ((-4 + 15 log 2) - ^ (557 + 840 log 2)p^ - ^ (2789 - 9660 log 2)pt + 0(pl)) 



'-"2) , 



A6,o(p2) = ^^^^^(2(47^661og2)p^-4(l-31og2)^3(89-1011og2)p^ + O(p^ 

A8.0(P2) = -^^^(4 + 2(l-141og2)p^_(i37-2241og2)p^ + 0(p^)), 

A9,o(p2) = -^(8(ll + 31og2) + 12(7-81og2)p^-3(91-2261og2)p^ + 0(p^)), 
7687rG V / 

Ak,,.2 (p2) = ^ ( - 12(1 - 2 log 2)pi + 3(25 - 36 log 2)pt + 0{pt)) , 

Aio,.a(p2) = ^(8(5-121og2)pi-3(91-1441og2)pt + C»(pi)), 

Ai„,;,2 (p2) = ^ ( - 90 + 4(29 + 60 log 2)p^ + (617 - 1080 log 2)p^ + 0(p^)) , 

A12..A (P2) = J ^ (20(-5 + 12 log 2)p2 + 5(71 - 192 log 2)p^ - (2713 - 3930 log 2)p^ + ©(pa)) 
5v37rG ^ ^ 

A12 x2 (p2) = ^= — f - 56(29 + 60 log 2)p2 + 840(1 + 16 log 2)p2 + (37403 - 51660 log 2)p2 + 0{p. 

35v37rG ^ 



and for the anomalous coefficients 

A2,a(P2) = ^^^(-24(l-log2) + 4(26-451og2)pi-(527-8221og2)pt + 0(pi: 

A5.4P2) = ^^ ( - 8(1 - 2 log 2)p2 + 16(5 - 4 log 2)pi - (389 - 566 log 2)p^ + 0(p5)) , 

A5.a(P2) = -^^ ( - 24(1 + 2 log 2)p2 - 8(17 - 48 log 2)p^ + 3(379 - 550 log 2)p^ + 0{pl 

A6.4P2) = ^^(-8(log2)p2-12(2-51og2)p^ + (174-2571og2)p^ + 0(p^)), 

A6.a(P2) = 7^ (24(5 + 2 log 2)p2 - 4(23 + 90 log 2)p-^ - (991 - 1494 log 2)p^ + 0{pl)) , 
487rG V / 

As,a(P2) = ^^(-241og2 + 4(25 + 421og2)pi + (199-7501og2)pt + C»(pi)). 

Note that in some cases the order 0(^2) vanishes, so that the corresponding expressions 
are valid up to Olpl)- 



C.25) 
C.26) 

C.27) 

C.28) 

C.29) 

C.30) 

C.31) 

C.32) 



C.33) 

C.34) 
C.35) 
C.36) 
C.37) 
C.38) 
C.39) 
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